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PHYSICAL SCIENCES 


EXISTENCE SEEN OF CYCLONIC EDDY IN NORTHERN EAST CHINA SEA 


Beijing HAIYANG YU HUZAO [OCEANOLOGIA ET LIMNOLOGIA SINICA] in Chinese Vol 16 
No 3, May 1985 pp 240-248 


[Article by Zaho Baoren [6392 0202 0088] of Institute of Oceanography at Qingdao, 
Chinese Academy of Sciences: "Recent Observation Evidence on the Presence of 
Cyclonic Eddy in the Northern Part of East China Sea"™] 


[Text] Abstract: This paper briefly summarizes the results of the surveys con- 
ducted in the cyclonic eddy zone in northern East China Sea in recent years. 
Based on the temperature field, mass field and flow field measured in these 
surveys, cyclonic eddy did exist in the area southwest of Cheju Island during 
the period. The residual current vector obtained in the survey was plotted 
with those measured in the past in the same season, resulting in a complete 
cyclonic eddy flow field. In addition, fundamental characteristics of the 
Huanghai Sea coastal current in northern East China Sea is also discussed. Our 
results further confirm that the diagnosis that this cyclonic eddy exists year 
round and its position is stable is correct. 


Based on temperature, salinity and current data gathered in July-August 1972, 

Hu Dunxin [5170 2415 2946] et al [2,3] pointed out in 1980 that a cyclonic eddy, 
approximately 60-90 nautical miles in diameter and 40 m in verticle scale, exists 
southwest of Cheju Island. In 1983, Mao Hanli [3029 3352 4409] et al further 
confirmed that this cyclonic eddy in northern East China Sea exists year round. 
In addition, the variation of this eddy in the summer season over the years was 
discussed based on the temperature field[5,6], All authors mentioned above be- 
lieves that the principal reason for the presence of this eddy is due to the 
horizontal velocity shear of the Huanghai Sea warm current and Huanghai Sea 
coastal current. 


Based on historic data, considerable information on the current north, east and 
south of the eddy is available. According to the vector diagram thus prepared, 
the flow field shows a tendency to rotate as a cyclonic eddy. Furthermore, 
these current vectors are consistent with other currents in the area such as the 
Taiwan Warm Current, Tsushima Warm Current and Huanghai Sea Warm Current. On 











the west side, i.e. the area reached by Hunghai Sea *cold Current and Coastal 
Current, however, there is relatively little information on the currents in 
that area. Hence, there are questions and doubts concerning whether there is 
a cyclonic eddy southwest of Cheju Island. To this end, we conducted three 
surveys on the hydrologic characteristics of the eddy in recent years. The 
focus of the measurements was to clarify the currents on the west side of the 
eddy. 
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Figure 1. The Survey Area!’ 
Key: 
® buoy; current measuring station in July 1982 
@ current measuring station in April 1934 
& buoy installed in joint Chinese-American survey in June 1980 
— cross-section position 


1) A buoy station installed during the joint Chinese-American survey in June 
1980 is also included in the fugure. 


The hydrologic cross-sections and current measuring statiors in these three 
surveys are shown in Figure 1. Based on an analysis, the flow of the Huanghai 


Sea Cold Current in the northern part of East China Sea was more or less under- 


stood. Once again, we proved that there is a cyclonic eddy in northern East 


China Sea southwest of Cheju Island. These observation results will be briefly 


discussed in this paper. 





*Survey Report No 1154, Institute of Oceanography, Chinese Academy of Sciences. 


This survey was performed under the guidance of Professor Mao Hanli himself 
and with the concern of Professor Guan Bingxian [4619 4426 6343]. Many com- 
rades from the technical survey group, hydrolic group and ocean current group 
participated in these three surveys. Figures in this paper were drawn by 
Comrade Du Meishan [2629 3270 1472]. The author wishes to express his grati- 
tude to all people mentioned above. 


Manuscript received on 26 January 1984. 











A. The Three Marine Surveys 
l. Results of Observation in July 1980 


The survey results centered around buoy arrays deployed in July 1980 were 
analyzed in detail in references [1,4]. References [5,6] also used some of the 
results. Hence, only a brief introduction is given here. Figures 2a and 2b 
show the temperature distribution and dynamic height distribution of the 10 m 
layer as measured in July 1980, respectively, as well as the average flow 
velocity vector of the 10 m layer as determined by an anchored array of buoys. 
In addition, we also included the average velocity vectors of the 14 m and 17 m 
layers obtained in the Chinese-American joint survey of the mouth of 

Changjiang River (from early to late June 1980)[7]. Based on Figure 2 we can 
see that (1) there is a cold center southwest of Cheju Island (125°30'E, 31°N). 
The lowest temperature is 19.2°C which is at least 2.7°C lower than those 
measured in surrounding stations. It was also found that the salinity in this 
cold center is 3.32 percent which is .15 percent higher than those measured in 
nearby stations. In the dynamic height plot, there is a corresponding closed 
contour. The Huanghai Sea Warm Current flows by the east and north sides of 
this contour, as shown in Figure 2 in reference [3], where there is an even 
higher dynamic height line. Hence, if we set up hydrological stations in this 
region, there will be more enclosed contours. The geostrophic cold water cir- 
culation will be clearer than that shown in Figure 2. (2) Based on the measured 
flow velocity vectors, the measured flow is very close to the geostrophic flow 
not only in direction but also in quantity. This fact shows that a geostrophic 
flow field can be used to approximately express the seawater flow in this 
region, at least in the summer. (3) As compared to the temperature chart 
(Figure 2a in this paper and Figure in reference [4]), the direction of the 
geostrophic current in the middle upper level and that of the measured residual 
current are in agreement with that of the isothermal lines in the area west of 
126°E. There is an overall tendency for the cold current of the Huanghai Sea to 
nove southward. This shows that there is indeed a cyclonic cold eddy in the 
northern part of the East China Sea southeast of Cheji Island. 
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Figure 2. Temperature, Dynamic Height Distribution and Measured Flow Speed at 

10 m Depth in July 1980. a. Temperature distribution; b. dynamic height dis- 
tribution (vectors shows are measured residual currents and values shown near the 
vectors are water depth levels.) 











2. Results ot Observation in Julv 1982? 


In order to further understand the hydrological structure of this eddy region 
in the summer and the flow west of this area, the survey vessel “Science 1" con- 
ducted a hydrological survey in Julv 1980. Figures 3a and b are the salinity 
distributions at 20 m and 35 m levels. Based on Figure 3, on the east side of 
the area surveyed, a cold center with an enclosed isotherm was found at 20 m. 
Its temperature is at least 5°C lower and its salinity is at least .05 percent 
higher than its surrounding environment. The center of this cold water region 
is located near 31°30'N and 125°45'E., At 35 m depth, this cold center has al- 
ready shifted northward by about 1 degree latitude to 32°30'N and 126°O0O'E. In 
addition, we can also see that there is a warm water tongue north of the cold 
center extending southwestward. This type of warm water tongue has been found 
several times before such as in July 1973 and June 19791). We believe that 
this is due to the branching of the Huanghai Warm Current before entering the 
Huanghai Sea. It is called the west branch of the Huanghai Warm Current. It 
often partially or completely cuts off the southbound Huanghai Cold Water to 
create an isolated cold water zone southwest of Cheju Island. 


Based on the residual current vectors shown in Figure 3 (Figure 3b shows vec- 
tors at 30 m depth), the flow direction is south on the west side of the cold 
center, which more or less agrees with the temperature distribution. Based 

on the data shown in Figure 3a on the flow on either side of the warm water 
tongue, the residual currents at 20 m and 30 m are mostly southwest, approxi- 
mately agreeing with the axis of the warm water tongue. Furthermore, the resi- 
dual currents measured at the two most north stations showed eastward flow. 
Based on the temperature field, this was a sign that the cold water eddy of the 
Huanghai Sea reached there. 





The cross-sectional diagram (Figure 4) shows that the high salinity cold water 
at the bottom is rising between 125° and 126°30"E. The cold water only reached 
20 m level, lower than previous years. 


In summary, in July 1982 the cross-sectional diagram clearly shows that cold 
water is rising. The horizontal chart shows a relatively strong high salinity 
cold center with an enclosed isotherm. Moreover, a warm water tongue extending 
southwestward was found on the northwest side of the cold center. Southward 
residual currents were measured west of the cold center. In addition, the 
direction of the residual current is more or less in agreement with the tempera- 
ture field. Based on these facts, the cyclonic eddy did exist southwest of 
Chenju during this survey. 


3. Results of Observation in July 1983 
Based on the data gathered in the spring of 1973 (April-May), Mao Hanlil5s6] 
et all believed that a cyclonic eddy southwest of Cheju Islands also existed 


in the spring. As shown in Figure 3 of references [5,6], the residual current 


1) Zhao Baoren, Current Structure in Cold Eddy Region in Northern East China 
Sea, [OCEAN AND LAKES], to be published. 








distribution measured in April-May 1973 is loca*ed south, east and north of the 
eddy. There is no data west of the region. For this reason, we conducted 
another hydrological survey in this area in April 1983. Five ocean current 

observation stations were set up on the west side in order to gather data re- 


garding the expansion of the Huanghai Cold Current into northern East China Sea ~ 
in the spring. 
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Figure 3. Horizontal Distribution of Temperature, Salinity and Residual Current 
in July 1982 
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Figure 4. Temperature Cross-section at 32°N in July 1982 
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Figure 5. Temperature and Salinity Distribution and Residuai Current in April 
1983 (same as in Figure 3) 


a. 20m level 
b. 35 m level 
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Figure 6. 32°N Cross-section Diagrams Obtained in April 1983 


a. temperature and salinity 
b. goestrophic flow speed 


Figures 5a and b show the distributions of temperature, salinity and residual 
current vector obtained in this survey at 20 m and 30 m. Based on the figure, 
there is a low temperature, low salinity tongue in the middle of this area 
extending from NNW toward SSE. This shows the tendency for the low salinity 





water of the Huanghai Sea to extend into northern East China Sea. On 
+} tA hj ~ :, hiol aitnitw . . 
elther side of this tongue, two lgnh Cemperature, hign sailnity tongues are 
Ye lew , a he ; T} 7 ~ ; Ass > + . 7 ~ + r 
extending northward. The west one is due to the northward extension of the 


Taiwan Warm Current and the east one is due to the northwestward extension of 


the Huanghai Warm Current. From Figure 5 we can see that the two high tempera- 
ture, high salinity tongues are approaching the cold water tongue north of 
> a T 


32°N. It is more apparent for the Huanghai Warm Current. Due to this effect, 
the low temperature, low salinity tongue south of 32°N appears as the head of 
a snake. Although temperature is rising in April, however, lamination effect 
is not obvious. From surface to bottom, the overall trend for horizontal dis- 
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The residual current vector obtained in this survey is relatively small, usually 
not more than 10 cm/s. Its direction is more or less in agreement with the 
tongue axis of temperature and salinity (see Figure 5). Dynamic calculation 
shows that the trend for geostrophic current in the cold water region agrees 
with that of the measured residuai current. Based on Figure 6a, a low tempera- 
ture, low salinity center exists at 125°-126°E. This is the bulk of the south- 
ward Huanghai Cold Current. From Figure 6b, the geostrophic current calculated 
to the right of this cold center moves southward and reaches a maximum. The 
calculated geostrophi< urrents on both sides are pointing northward. Thus, 


the geostrophic current pattern agrees with that based on the temperature and 


salinity tongue distribution. Similar tendency was found along several cross- 
sections in the south. The difference is that the center of the southward 
geostrophic current also moves southeastward as the cold center moves in the 


same direction. 
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Figure 7 shows the dynamic height obtained in this survey at 20 m and 30 a, 
Due to lack of data, only a partial elliptical cyclonic circulation is drawn 
southwest of Cheju. Considering the flow of the Huanghai Warm Current south- 
west of Cheju in the spring, it is possible to complete the ellipse. Hence, 
during this survey period the cyclonic eddy also existed. 


II. Huanghai Coastal Current in Northern East China Sea 


There is a wide range of low salinity cold water in northern East China Sea. 
As for the important questions such as the dynamic nature of the area, how 
does it reach northern East China Sea, and how fast does it flow, some funda- 
mental understanding has been obtained based on analyses of the three surveys. 


Based on Figure 2, the summer Huanghai Coastal Current flows southeastward with 
the cold water density circulation in southern Huanghai Sea. When it reaches 
the 32°N line east of 124°E, it is deflected toward the southeast on the west 
side of the cyclonic eddy described above. Then, it merges with the Taiwan 
Warm Current south of 31°N and east of 124°E and then moves toward southeast. 
Based on residual current vector measured, Huanghai Coastal Current reaches 
10-15 m/s after it gets to East China Sea. It mixes with the surrounding sea- 
water continuously. First, it mixes with the warm water carried by the western 
branch of the Huanghai Warm current near 32°N and then merges with Taiwan warm 
water further south (south of 31°N). In addition, the upper layer is mixed with 
the fresh water of the Changjiang River. 


As described earlier, cold water is delivered to East China Sea by the spring 
Huanghai Coastal Current in the direction of the tongue axis (southeast). At 
this time, the direction of Huanghai Coastal Current is monotonous. It fiows 
consistently southeastward from southern Huanghai to northern East China Sea, 
Based on the data and dynamic calculation, its flow speed can also reach 

10 cm/s in the spring. 


Based on a comparison of dynamic calculation and residual current measurement, 
in spring and summer, the Huanghai Cuvastal Current is basically geostrophically 
balanced as it moves south. But, the situation in spring is different from 

that in the summer. In spring, the temperature difference between the south- 
bound cold water and the surrounding water i.. relatively small. Because of its 
apparent low salinity, although the Huanghai) Coastal Current is geostrophic, it 
moves along the axis of the temperature and salinity tongue. In summer, it 
moves southward along the west side of the cold center. At this time, the 
isotherms approximately reflects the geostrophic flow. We do not have any 
residual current data in fall and winter seasons. However, based on the geo- 
strophic flow, its distribution in winter is more or less similar to those shown 
in Figures 6 and 7. In the fall, before the end of October, it is approximately 
the came as that in the summer. Thus, the coastal current in winter is expected 
to be the same as spring and the situation in the fall is similar to that in 
swamer. 


III. Flow Field Measured in Eddy Region 


Another important question to be explained is that whether the flow velocity 
vectors measured on the west side of the eddy agrees with other vectors observed 








during the same period in the past. Figures 8 and 9 show the results combin- 
ing the re:.dual current vectors observed in the summer (solid lines) and 

those obtained in July-August 1972 (dotted lines). Figures 8a and 9a show the 
residual current distribution at 10 m and Figures 8b and 9b: w that at 20 m. 
Based on these four figures, the residual current vectors ob ‘d in the summer 
of 1930 and 1982, together with those measured nistorically in . 2 same period, 
form a complete cyclonic eddy motion. The center of the eddy is approximately 
located at 126°E and 31°30'N. The diameter of the eddy is 2 lattitude units. 

It should be pointed out that this situation not only occurs at 10 m and 20 mn. 
From 5 m to 50 m, similar trend exists. 
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Figure 8. Residual Current Together with Data Obtained in Summer 1972 


a. 10m layer 
b. 20m layer 
(Dotted vectors are 1972 summer data.) 
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Figure 9, Measured Residual Current Combined with Summer 1973 Data. 


a. 10m layer 9 
b. 20m layer 


(Dotted lines represent 1973 
summer data.) 
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Figure 10. Measured Residual Current Combined with Spring 1973 Data 


a. 20m layer 
b. 30m layer 
(Dotted lines represent spring 1973 data.) 


Figure 10 shows the combination of the residual currents measured in April 1983 
and April-May 1973. Based on the figure, similar to the summer flow field, the 
combination of flow velocities measured over the years in the spring form a 
complete cyclonic flow field. The difference is that the spring flow field 
appears as an ellipse in the north-south direction which is consistent with 
the distribution curve shown in Figure 7. Furthermore, the center of the eddy 
is shifted southeastward by a half latitude unit. The spring eddy thickness is 
over 50 m. 


IV. Conclusions 


Major findings of hydrological surveys conducted in northern East China Sea 
have been described above. The flow of Huanghai Coastal Current into northern 
East China Sea has been explained. We also pointed out that a cyclonic eddy 
southwest of Cheju Island was present during every survey. By plotting all the 
historic residual current data together, the result shows a complete cyclonic 
eddy in that region. This cyclonic eddy motion indeed exists in the region. 
eerie ae | we proved that the oe cm) ees reached by Hu Dunxin [7579 2414 
2946] et al 3] and Mao Hanli et al [»®! that a semi-permanent year round 
cyclonic eddy exists southwest of Cheju Island is correct. 
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PHYSICAL SCIENCES 


SPECTRAL ANALYSIS OF WAVES IN GULF OF TONKIN DESCRIBED 


Beijing NANHAI HAIYANG KEXUE JIKAN [NANHAI STUDIA MARINA SINICA] in Chinese 
No 2, Oct 81 pp 79-92 


[Article by Wang Wenzhi [3769 2429 6347], Nanhai Oceanography Institute, 
Chinese Academy of Sciences, "A Preliminary Analysis of the Wave Spectrum in 
the Gulf of Tonkin during the Passing of a Cold Wave"; paper received 

10 October 1978] 


[Text] Although lots have been accomplished in the study of sea-wave spectrum, 
its spectrum form has not yet been unified. The various spectrum representa- 
tions that have been proposed give the results that are quite different from 
one another. Therefore, it is extremely important to analyze the observed wave 
spectra of a selected area of the sea. On the one hand, it will produce know- 
ledge about the composition of the wave spectrum of a certain sea zone, which 
will be useful for areas such as constructions at sea, ship-building and wave 
forecast. On the other, it will provide information helpful to the comprehen- 
sive efforts in the theoretical study of sea-wave spectrum, 


Spectrum calculations were carried out on 12 wave records from a fixed site in 
the Gulf during the passing of a stormy cold wave and a preliminary analysis and 
discussions of the results are presented in this article. Among the 12 calcu- 
lated correlation function graphs, some are not perfect. Either they don't 
oscillate along the time axis or they are far from approaching zero. These are 
probably due to insufficient measuring time and the presence of the "zero-line 
drifting" phenomenon, which will be discussed in this article. 


1. Source of Data 


The wave data used in this article were measured by a wave gauge on board a 
Chinese ship at a hydrological station (109°15'E and 20°45'N) in the Gulf of 
Tonkin between 23:30 of 16 January 1973 and 14:00 of 17 January 1973. For 
practical reasons, each measurement lasted only about 10 minutes. At the same 
time waves were measured, wind speed, wind direction, and water level were 
also recorded. Wind speed was recorded by an automatic anemometer, which re- 
cords the average wind speed at 2 minute intervals. Table 1 lists the wind 
speed, wind direction and water level changes recorded during the passing of 
this storm and Figure 1 shows the changes of wind speed and tide level. 


12 





Figure 1. Changes of wind speed and water level at a hydrological station in 
the Gulf of Tonkin during a storm on 17 January 1973. 
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Table 1. Recorded wind and water level data 


























q ani » A, 8 cR q ARES AKAP SKE MK 
(8 ,i:5) (%/®) (28) e (OR) Fo Cm) 
16, 23:30 12.0 N 77 20 
16, 23:52 20 
17, 92:00 11.0 NNE 95 20 26.7 
17 02:00 13.5 NNE 95 20 27 .0 
17,04:00 15.6 NNE 95 20 26.8 
17 05:00 13.5 NE 20 27 .3 
17, 06:00 14.0 NNE 95 20 27.3 
17, 08:00 16.5 NNW 20 25.8 
17 09:00 20 
17, 11:00 15.0 N ? 20 
17, 12:00 14.0 N 7 20 23.0 
17, 14:00 13.0 N 77 20 23.$ 
Key: 

a. Observation Time (day, hr:min) d. Fetch length (Km) 

b. Wind speed (m/s) e. Average depth within the fetch (m) 

c. Wind direction f. Measured water level (m) 


(Subsequently, the data will be designated by their observation times.) 

2. Methods and Results of the Spectrum Calculations 

We treated the sequences of wave observed as a quasi-stable process. The con- 
tinuous wave records were scattered by a time interval of at = 0.375 s to get a 
series of samples from the time sequence $xi=x(iat), i = 1,2,**+, N& . The wave 


spectra were then calculated in four stages, 


(i) Calculations of the average values 


S* x, (1) 
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(2) Calculations of the correlation functions 


C,=C(ra:)= ~—> (2, — 224, — 7) 


iz|i 








g LS by K*x? 
z5 — 
N-rm@ 
y™=0,1,2,°**,™ (2) 


> 


where K is a factor for the conversion of height on chart paper into actual 
height. 


(3) Calculations of the unrefined wave spectrum 


a- i 
V(a) = |Co+ 2S) C,cos Sh + Cacos x4] 


Tai 


in which wo, = hie = ——h 
mit 


hm0,1,2,°¢*,m (3) 


(4) Calculations of the refined wave spectrum 


S€ey) = 0.547 Cea) + 0.4600 @,) 
SC ao, Y = 0.94F (co, ) + 0.238 PC co.-1) + PCs, )) 
hmm l,2,°++,m—1 


SC.) = 0.54V(w.) + 0.46V(w.-,) (4) 


Actual calculations were carried out on a model DJS-21 computer. The data 

for each time sequence are shown in Table 2. Through the calculations, 12 
correlation functions and spectra were obtained, 4 of them are shown together 
in Figure 2. To make it easier for comparison, we divided the 12 spectra 

into 2 groups: those before 170800 in one group and those after in the other. 
They are shown as a group in Figures 3 and 4. 


In the program for wave spectrum calculations, we also included the calcula- 
tions of three relevent features: 


(1) Calculations of the zero-, 2nd-, and 4th-order moments of the spectrum, 
Mo, M2, M4, and the spectrum width coefficient &. The equations are as 
follows: 


Mom (=) (1 + $C coy) + 4+ S(em) + 2+ S(coq) + +++ 
edge S( ate, ) + 1+ SCa,)) 


Me () OF + SCem) + 4+ 1+ SC) 


+ 2+ 2+ SCco) + 22) 4+ (m— LP + Sides) 
+1: m+ S(@,)) 
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1/ «= \%, ¢.e ‘ 
\ [1+ O*- Sag) + 4° 1S - Sim) 


+2-2 * Senn) + -+4-(m—1)- S(o,-1) 
(5) 





(6) 


Wan 
and are listed in Table 3. The calculated 


—E= — 


The calculated values of Mo,M2,M4, 
Mo should equal to 2Co. 














“0.24+ 


Figure 2. The correlation functions of four sets of reccrds. 








Table 2. Related parameters in the spectral analysis of the recorded data. 











@ Syo |5os (®) ye) | ae RE m 4 
162330 0.373 1124 7.0 42 3 
162352 0.375 1440 9.0 $6 $1 
170000 0.375 1468 9.2 145 56 $1 
170200 0.375 2048 12.8 153 84 “8 
170400 0.375 2048 12.8 134 84 “8 
170500 0.375 2048 12.8 150 84 48 
170600 0.373 1440 9.0 100 56 s1 
170800 0.373 2048 12.8 134 8 8 
170900 0.375 1300 9.4 109 56 $3 
171100 0.375 2048 12.8 147 84 48 
171200 0.375 1312 9.5 110 $6 $4 
171499 0.375 2048 12.8 151 4 8 























t--time interval; N--sample capacity of a time sequence; T = N- t--total dura- 
tion of a time sequence; m--# of data point for correlation function and spec- 
trum; df = 2[(N/m-1/4]--degree of freedom of the calculated values. Values of 
50-60 are appropriate for wave spectrum analysis. 


Table 3. Calculated values of moment and spectrum width coefficient. 








ae 8 Mm | M, (m= S* u,(m> S-#) ; 
No. Ce") (m2) 8") RO~) 
162330 0.026 0.240 7.408 0.834 
162382 0.106 0.356 4.014 0.837 
° 170000 0.392 0.430 3.490 0.954 
170200 0.278 0.578 3.852 0.829 
170400 0.430 0.940 9.826 0.912 
170500 0.328 0.688 4.946 0.841 
170600 0.512 0.888 4.788 0.823 
170800 0.702 1.240 9.964 0.883 
170900 0.588 0.972 5.394 0.838 
171100 0.292 0.578 4.862 0.875 
171200 0.428 0.854 8.750 0.897 
171400 0. 360 0.870 11.006 0.900 

















(2) Examing the changes of the correlation functions. Given§> 0 and a posi- 
tive integer dm, a number m is selected by computer such that the absolute 
values of the finald m correlation function values Cj do not exceed d , that 
is 


ic,| <2 py=em—8,+1,m— 8, + 2,°*+,m 
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Figure 3. The Spectra of 162330-170800. 
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Take 170800 as an example. When = 0.1 andj m = 20, the result is m = 38; 
Wheng = 0.01 and dm = 20, it is m> 222. It is clear that even with a very 
large m the oscillation amplitude of the correlation functions is not completely 
reaching zero. 


(3) Dividing each time sequence into four sections and calculating their 
averages. Given the averages of the four sections, comparisons of them with 
the total average of the time sequence can be made. We define the maximum de- 
viation of the section averages from the total average as 

max{ s;— 2} 
i. 


i= {—1,2,3,4 (7) 





z 


From the values of the maximum deviation we can determine whether the average 
values of a time sequence are constant or near constant, which would lend 
support to our treating the sequences of wave as a stable process. Of course, 
this is not the only criteriun. To be a stable process, the condition that 
“correlation functiou is related only with time interval" has to be satisfied 
also. We did not check if this condition was met. 


The 4 values for the 12 time sequences are listed in Table 4. Except for 
170000, the A values of all sequences lie between 0.8-6.7 percent. The devia- 
tions are very small. This provides some ground in our treating approximately 
the sequences of wave as a stable process. 


Table 4, Maximum deviations of section averages from the total average. 

















@ No. 8 & (%) @ No. & | AR (%) 
162330 2.1 170600 $.4 
162352 3. | 170800 +4 
170000 31.7 170900 3.9 
170200 4.0 171100 3.2 
170400 6.7 171200 $$ 
170500 0.8 ’ 171400 1.4 





3. Analysis and Discussions of the Results 


From their appearances, our graphs show only the energy distributions of the 
component waves within a frequency range and particularly that their peaks 
all fall within a narrow frequency range. It can be seen from Figure 3 that 
the curves become steeper as wind speed, wind duration and fetch increase and 
the total area under the curve also increases. Also the spectral peaks span 
a wider frequency range and the corresponding wave height and period of the 
wind-wave increase accordingly. Meanwhile, during the over 2 hours period 
between 162330 and 170200, the peaks shift rapidly toward low frequency end 
so the spectra look quite different at low frequency region but very similar 
at high frequency region. But the phenomenon becomes insignificant after 
170200. Though the peaks still shift toward low frequency end, it happens at 
a far slower speed. Not only the curves are very similar, but the maxima of 
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the spectra occur at nearly the same frequency. The values of later-stage 
spectra are generally greater, at all frequencies, than the corresponding ones 
of early-stage spectra. It can be seen from Figure 4 that the peaks diminish 





ind the areas under them decrease as wind speed and fetch decrease. But the 
frequency range of these peaks changes very little and the decrease of the 
-orresponding wave height is more noticeable than that of the period. 


There is a common feature in this two groups of spectral graphs, namely the 
spectrum values do not approach zero at the very low frequency end of these 
spectra (nearw = 0). This is physically impossible. To certain extent, this 
phenomenon can be attributed to calculation errors. And these errors have to 
do with the errors introduced in our determination of the zero-line of a wave 
profile. In other words, the time sequences are not completely centered. In 
the spectral analysis of these wave data, we simply used the total averages of 
the time sequences and that is not good enough. Due to the effect of tidal 
fluctuations, the effect of “zero-line drifting” in the recorded data, or 
other errors in calculation, the time sequences are not completely centered. 
The hint and extent that the 12 sets of the data analyzed are not completely 
‘entered can be seen from theirA values, correlation function graphs and 
spectra. Three key data are listed in Table 5. It reveals that the spectrum 
value atW= 0 (i.e. the nondimensional quantity S(0)MoT in table 5) is not 
issociated with wave strength but is dependent on the A value and the position 
of the correlation function curve. To be specific, the greater the and the 
farther the correlation function curve is away from the time axis, the greater 
the spectrum value at Ww = 0 and the worse it is off the center. 

Table 5. The effect of off-center on the correlation function graph and 
spectrum, 
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When we select the scattered dat 
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2d cata S reco we assuMec 1 
zero-line of the measured curve wa Straight lime and the readings were made 
based on a line parallel to the assumed zero-line (such as baseline I in Figur: 
5). When the zero-line of the me d rve drifts (line the curve in Figure 
5), then it is unlikely to be completely centered. If we re-examine the cir- 
cumstances under which the curves were measured and draw as closely as, possibls« 
the actual zero-line (not necessarily a straight line), and then select a 
baseline II that is parallel to it (as shown in Figure 5) for data readings, 
the extent that the time sequences are centered can be greatly improved. We 
purposefully applied the method / nd 1/0600, which have the largest 
zero-line drifting, to modify the data and carried out similar calculations. 
The results are shown in Figure 6 and Table 6. When compared with the graphs 
of 170000 and 170400 in Figures 2, 3 and Table 5, it is obvious that there is 
a great improvement. In the part that w/0.5, the spectral density graphs al- 
most completely coincide with each other. Also for the part that «<0.5, par- 
ticularly near where & =0, the spectral density values are very close to zero. 
The same for the correiation fucntions. heir graphs agree totally with each 
other and they oscillate around the time axis and are more symmetrical. In 
short, they are closer to reality. 

ture 5, Schemat i plot of baseline I an li, 
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Table 6. The improvements in correlation function graphs and spectra after the 
degree of centering is improved. 
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Figure 6-a. The calculated correlation 
using the data based on baseline II. 
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Figure 6<b. The calculated spectra of 170000 and 170400 using the data based 
on baseline II. 
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in order to accurately estimate the total area of the spectral peaks (equiva- 
lent to the total energy of the wind-wave), we extracted AMo, the area corres- 
ponding to the region of 0< uw<0.25 s~l, from the total area Mo to get the modi- 
fied total area E of the spectral peaks. The results are listed in table 7, 
Figure 7 depicts the changes of the total area based on the data in table 7. 


Table 7. The calcualted areas of the spectral peaks. 











’ VP 4 

“ee ] ese dey | eNOS os ES ces) 
162330 0.026 0.0042 0.022 170600 0.512 0.0398 0.472 
162352 0.106 0.0224 0.084 | 170800 0.702 0.0040 0.698 
170000 0.592 0.4734 0.119 | 170900 0.588 0.0398 0.548 
170200 0.278 0.0146 0.263 171100 0.292 0.0202 0.272 
170400 0.430 0.0738 0.356 || 171200 0.428 0.0308 0.397 
170500 0.328 0.0080 0.320 171400 0.360 | 0.0102 0.350 























Figure 7. Changes in spectral peak areas, 
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We have mentioned earlier that, as shown in Figure 3, the frequency range that 
covers the peak and major energy distributions of the wind-wave spectrum 

shifts from high toward low frequency end. But after 170200, the change be- 
comes less obvious. The major distributions of energy are concentrated within 
the frequency range of [0.65, 1.75]. Table 8 lists the exact values of 

[S@v) ]max of the individual spectrum before 170800, which reach the maximum 
near the frequency Wm. The frequency ranges that cover the major distributions 
of energy are also listed. 
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The last column in Table 3 were calculated from eq.(6). This is one way to ob- 
tain the spectrum width coefficient € from the internal structures of the wave. 
One can also get € from the shape of the measured curves for comparison. 


Table 8. The distribution characteristics of the wind-wave energy with respect 
to frequency 


























*.% *'D) or) - 
(m2.S (s7+) ASHES") ci S4ERRHESKd 
162330 0.018 2.20 (0.8,5.8) 8s 
162352 0.120 1.70 (1.35, 2.55) 86 
170000 0.190 1.65 (1.35,2.10) 8s 
170200 0.468 1.05 (0.65,1.75) 8s 
170400 0.676 1.01 (0.65,1.75) 8s 
170500 0.664 1.06 (0.65,1.75) 8s 
170600 0.796 1.11 (0.65, 1.75) 85 
170800 1.300 1.07 (0.65,1.75) 85 
Key: 
a. No 


b. Frequency range that covers the major distribution of energy 
c. Frequency range (s~1l) 
d. Percentage of the total energy 


If one defines the average period of the measured curve x(t) as Te = t'/Ne, 
and the average period of dx(t)/dt as Tc = T'/Nc where t' is the time duration 
of x(t), Ne is the number of times that x(t) reaches zero within time t' and 


Nc is the number of peaks (the maxima of x(t) values) within the t', it can 
be proved that Tc/T,)2 = M_2/MoM4. Therefore, eq.(6) can be rewritten as 


$B 
en fi-(%) (8) 


The € values listed in table 9 were calculated from eq.(8). 





or 





It is known that ¢ is a distinctive quantity depicting the degree the wind-wave 
energy is concentrated with respect to frequency and that 02€<1. When € is 
closer to 0, the frequency range in which the spectrum densities converge be- 
comes narrower. In the literature [3], those with ¢< 0.4 were referred to as 
narrow sepctrum and those with £? 0.4 wide spectrum, As evident from Tables 

3 and 9, ovr calculations show that the results from eqs.(6) and (8) are quite 
different, although these equations are theoretically identical. Generally, 
the calculations of moment involved in eq.(6) are inaccurate, particularly for 
higher order moments M2 and M4. Therefore, the calculation results of eq. (8) 
seem more reliable. 
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Table 9. The calculated values based on the shape feature of the measured 








curves. 
| No = Ny i N, 8 
162330 131 147 0.454 
162352 147 166 0.465 
170900 . 148 157 0.334 
170200 176 204 0.506 . 
170400 180 196 0.396 
170500 191 205 0.363 
170600 120 131 0.401 
170800 172 193 0.454 
170900 103 116 0.460 
171100 iSS 175 0.464 
171200 119 136 0.484 
171400 167 190 0.477 














The spectrum is a function of such wave elements as wave height and average 
period. For example, 


Amwi77/8, Ayam 283VE 


A, = 3.60 /E (9) 
?=2z Ms (10) 
M, 


Table 10 gives the comparison of the measured values of wave height and period 
with those calculated from eqs. (9) and (10), in which the “measured values" 
were obtained by statistical means from the measured curves. It can be seen 
that the measured and calculated wave heights agree well while the calculated 
periods are significantly lower than the observed, which is probably due to the 
larger errors in the second-order moment M2. 


Table 10. Comparison of the calculated and observed values of wave height and 














period. 
7 ve | BORA) BCR Huw (® t (p) G 
9 ceo | - 
wa LSP | AP Ladle [fe [ter re ne dy re 
762330 | 0.022 | 0.148 | 0.26 0.42 0.54 0.026 | 0.240 | 2.6 
162382 | 0.084 | 0.289 | 0.31 0.82 1.05 0.106 | 0.386 | 3.4 


170000 | 0.119 | 0.343 } 0.61 | 0.59 | 0.97 | 0.97 1.24 1.31 | 0.592 | 0.430 7.4 3.7 
170200 | 0.263 | 0.512 | 0.91 | 0.97 | 1.45 | 1.48 1.85 1.99 | 0.278 | 0.578 4.4 $.0 
170400 | 0.336 | 0.995 | 1.06 | 1.08 | 1.69) 1.65 2.15 | 2.18 | 0.430 | 0.840 4.3 3.0 
170800 | 0.320 | 0.564 | 1.00) 1.06 | 1.60) 1.64) 2.04 2.17 | 0.328 | 0.688 4.3 $.1 
170600 | 0.472 | 0.685 | 1.22 | 1.30) 1.94) 2.01 | 2.48 2.52 | 0.512 | 0.888 4.8 3.4 
170800 | 0.698 | 0.833 | 1.48 | 1.30 | 2.36) 2.30) 3.01 2.90 | 0.702 | 1.20 4.7 5.0 
170900 | 0.348 | 0.739 | 1.31 | 1.30] 2.10) 2.06 | 2.67 2.60 | 0.588 | 0.972 4.9 $.2 
171100 | 0.272 | 0.$20 | 0.92 | 0.98) 1.48) 1.57) 1.88 1.99 | 0.292 | 0.578 4.5 $.2 
171200 | 0.397 | 0.629 | 1.12 | 1.14) 1.76] 1.84) 2.27 2.31 | 0.428 | 0.854 4.4 $.2 
171400 | 0.350 | 0.589 | 1.05 | 1.04] 1.67] 1.65) 2.13 2.22 | 0.360 | 0.870 4.0 $.1 
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4. Conclusions 


The spectral analysis of 12 sets of observed wind-wave data during the passing 
of a storm was carried out and certain insights about the wind-wave spectrum 

of the Gulf of Tonkin during the passing of a cold wave were gained. However, 
some of the information needs to be verified with additional observations. As 
to the collected wave data, we have to further improve the aspect of data treat- 
ment and further refine our analysis methods. First, zero-line drifting exists 
to a varying extent in the observed wave data. In our attempts to eliminate 
this effect, we should try various approaches and compare their results; secondly, 
lots of personal judgment were involved in our selections of m values in eqs.(2), 
(3) and (4). Different m values will give different spectra. The greater the 

m, the sharper the spectrum peak. So which m value will give a corresponding 
spectrum that matches the real spectrum? Thirdly, regarding the problem that 
the two spectrum width coefficients as calculated fromthe internal structure and 
external characteristics of the wind-wave are different, how can they be unified 
on the basis of their formulas and the physical meaning of& ? Fourthly, the 
data collecting during the passing of the storm was rather thorough. We ought 

to investigate further the relationship of the spectrum change with wind speed, 
wind duration, etc. In summary, the analysis of these wave data has to be 
continued. 


Acknowledgement: After the manuscript was written in December 1973, Professors 
Wen Shengchang [2429 5110 1603] and Yu Zhouwen [0151 1352 2429! of Shandong 
College of Oceanography provided valuable advices and enthusiastic guidance. 
Recently, additions and revisions have been made. We want to express our sin- 
cere thanks to Prof Yu Zhouwen for reading the final manuscript and Comrades 
Chen Junghang [7115 0193 2490], Zhang Jinghan [1728 4843 3352], Lan Changhua 
(5663 2490 5478], Sui Shifeng [7131 0013 1496], Jia Shufen [6328 3219 5358) and 
Xian Huilan [3156 1920 5695] of the Wave Study Group of the Institute for their 
participation in data handling and advices. 


REFERENCES 
[1] Wen Shengchang, "Wave Theory," Shangdong People Press, 1964. 


[2] Blackman, R.B. and J.W. Tukey, 1959, "The Measurement of Power Spectrum 
from the Point of View of Communications Engineering," New York, Dover. 


12922/1223? 
CSO: 4008/1052 


26 














PHYSICAL SCIENCES 
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Guangzhou; paper received 12 October 1982, finalized 4 April 1983] 


[Text] Abstract: The use of acoustical methods to obtain information on 
the physical-mechanical properties of the ocean floor wrgere 4c 3 | layer has 
been one of the fast developing ocean floor research methods'*"*’ in recent 
years. Observations during the last 10 to 20 years using electron micro- 
scope technology have led to the discovery of a few important microscopic 
characteristics of sedimentation; for example, the condition of granular 
sediment buildup, the pore structure, and aspect ratio of pores, etc. All 
these have effects which cannot be overlooked when considering the propaga- 
tion characteristics of sound in ocean floor sediments, particularly in 
terms of the velocity of sound. In order to clarify the structure of ocean 
floor sedimentation and the physical mechanism of the mechanical factors that 
influence the propagation characteristics of sound, the authors investigated 
the relationship between the microscopic structure of sediment and macro- 
scopic acoustical phenomena. Also, taking one representative area of the 
continental shelf in the northwest section of the South China Sea as an 
example, regression analysis was carried out on a set of test data from 17 
test sites. Letting the unidimensional compressive strength, qy, be equiva- 
lent to a factor which comprehensively expresses the macroscopic effect of 
the granular buildup and pore structure of ocean floor sediment, the 
porosity, n, was used to calculate a supplementary revised term for the 
speed of sound, c, in the sediment layer. The following empirical formula 
was obtained: 


f(c,m,q.) =C—4195 +90, 5833n—0, 7695n*° +94 69689, — 64 4603gi = 0, 
The formula raises the calculation precision of the often used formulae (for 
example, the Anderson formulal5]) which do not consider pore structure; and, 


at the same time, it eliminates the necessity of carrying out a statistical 
survey of corrected sound velocity calculations on microscopic pores. 
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1. Theory and Present Probiems 


Experimental resuits and theoretical analysis point out that, when sound 
waves passing through a type of fluid such as ocean fioor sediment become 
oversaturated in the medium contained in the pores, their propagation 
velocity should be influenced by the fluid portion and the solid portion of 
the medium. M.A. Biot, in a series of papers|6], proposed a theory for 
sound propagation in oversaturated pore-contained fluids. This theory was 
also the first to point out the second power relationship of the speed (i.e. 
sound velocity, c): 


c=(P+20+R)/p, (1) 
where P, Q and R are the Biot coefficients, namely: 


_ (a—1)* C,+(n—-1)C.,4+C a, 
P= C(C.—-C..)+C,.C., +4/3u, 
_ _ (n—:}\C,+C., 
Q=—" OC ~C.,)+C,C., * 


R= eC: <a 

(C.—C.,)+C,C., 

in which, C; is the volume compression coefficient of the total sedimentary 
framework; Csp is the solid granular compression coefficient; Cp is the pore 
compression coefficient; Cy is the fluid compression coefficient; n is the 
porosity; and Ct = nCp + Cgp. 


Althouch Biot's considerations may be all-encompassing in theory, in prac- 
tical application these coefficients have almost no possibilities for total 
measurement. Empirical results also indicate clearly that it is not neces- 
sary to measure all these coefficients. For this reason, the simple M.R. 
Wyllie model and the "Wyllie sound time averaging formula" for porosity in 
oil bearing sediment and sound velocity related calculations are generally 
adopted as a basis. 


An example of Wyllie's sound time calculation model is shown in Figure ]. 
It indicates the two substances (a solid consisting granular solids and a 
liquid or gaseous material consisting of fluid saturated with granular 
particles) composing the porous media. The propagation velocity of longi- 
tudinal waves is equivalent to the velocity through the two types of media 
according to their volumes when totally separated. Therefore, the sound 
time, t = l/c 1), in passing through the composite medium should equal the 
sound time of passing through the solid portion plus the sound time of 
passing through the fluid portion: 


t = (t¢ - tma)n + tma (2) 








1) 1 has a unit of length, dimension: L. 











in which, t¢ is the value of the sound time in the fluid within the sediment 
pores and tma is the value in the solid framework of the sediment. 








Figure 1. Fundamentals of the Wyllie sound time averaging formula model; 
the slanted line portion denotes the solid medium, the remainder 
is the fluid medium; arrowheads indicate the direction of sound 
wave propagation. 


From this, we obtain the relationship between the sediment porosity, n, and 
the sound velocity, c: 


! 


Ce alc.s—e,) +e, 


—(Ce.°C;), (3) 


where the subscripts of c are the same as in Equation (2). 


These concepts have found extensive application over the last 20 years in 
areas of petroleum geology and acoustical well logging; moreover, they have 
become the theoretical basis for applied formulae for sound velocity/porosity 
calculations often used in sediment acoustics (for example, the Anderson 
formulal5] and formulae quoted in other papers|2,8]), 


But, the Wyllie model and its resulting relationships also fail to take into 
account the natural condition of the pattern of buildup of granular sediment 
and the state of mechanical coupling. They especially do not consider the 
form of the sediment porosity and the influence of the dominant orientation 
of the porosity. Generally speaking, the form of the sediment porosity is a 
random distribution, but in view of all existing primary and secondary 
structures like microstrata, microcrevices, spaces between layers, or con- 
ditions of statistically dominant orientation of pores of long axis 

(Figure 2), there appears to be a relatively large discrepancy between the 
value of sound velocity as derived from the Wyllie formula and that from 
actual measurements. Biot's theoretical model and Equation (1) have problems 
similar to these. Because of this, in 1976, M. Nafi Toksoz et al. furthered 
the investigation into these problems. In calculations based on porous media 
of equivalent elasticity constant, they showed evidence of porous media of 
non-spherical pores with different aspect ratios. The velocity of longitudi- 
nal and transverse waves showed very great differences as compared to porous 
media of near spherical pores; among the differences, the effect of 








longitudinal waves were greater than that of the transverse waves|l9]. 
Because of this, they proposed the concept of a pore shape chart and also 
proposed the necessity of determining the aspect ratio of pores, a, to 

carry out corrections on sound velocity calculations. Wei Moan [7614 1075 
4152] and Zhu Zhengya [2612 2973 0008]'+°) furthered the advances by obtain- 
ing a pore structure factor, K, on the basis of the determination and 
analytical calculation of actual microscopic patterns of pores. This 
describes, in the longitudinal wave direction, the effect of pore shape and 
effective pore diameter on sound velocity. 














Figure 2. Structure of granular sediment buildup and inter-granular pore 
Orientation; arrowheads indicate the direction of sound wave 
propagation 
(a) ideal isotropic structure of sediment and particle axes; pores 
have no dominant orientation; sound velocity in each direction is 
equivalent 
(b) the long axes of the particles and the pore shapes have a 
definite dominant orientation; sound velocity varies in each 
direction 


Nevertheless, regardless of Toksoz's pore shape spectrum and the aspect ratio, 
or the pore structure factor, K, correction applied by Wei Moan and Zhu 
Zhengya to the Wyllie model, each must make accurate measurements of pore 
size and shape, again carry out statistical operations, and definitely 
include the differences of the aspect ratios. This requires the use of 

water or oil saturation of future samples, or use of vacuum manometer 
methodology, even to the point of using the method of vacuum impregnation of 
colored organic glasses, and the microscopic examination of pore shapes, and 
more calculations of their horizontal and vertical contrast, a, or pore 
structure parameter, K. This way of doing things is not only technologically 
tedious, it's also difficult to obtain solutions to practical problems. 
Moreover, one must guess at the porosity, the aspect ratio, and the pore 
structure in a statistical distribution which is homogeneous, particularly 

in regard to the effect of a closed pore medium which does not allow a way 

to carry out an observation. All this limits the practical use of the theory 
and also limits further advances in this aspect of the research. 
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If. Sediment Mechanical Parameters and Velocity of Sound Data 


Since 1980, at a 2100 square kilometer area of the continental shelf in the 
northwest section of the South China Sea (60 kilometers east-west by 35 
kilometers north-south), we carried out a series of gravitational column 
samplings of ocean floor sedimentation, included work on an ocean floor 
acoustics investigation, and researched the physical-mechanical characteristics 
of this ocean bed material. Figure 3 shows the distribution of the column 

sam ling points and the sound velocity measuring points within the investi- 
gated area of ocean. Analytical test results of each earth-mechanical 
parameter and elasticity parameter of this area of ocean floor sediment are 
listed in Table 1 and Table 2. 














Figure 3. Distribution of test points in the experimental ocean area 


The sediment mechanical parameters listed in Table 1 are all according to the 
People's Republic of China 1979 issue of SDS 01-79 soil mechanics regulations 
for carrying out testing and test conclusions; errors are all compatible with 
regulation requirements. The included specific gravity, G, is the specific 
gravity as measured by the bottle method; experimental results indicate that 
in the sediment in this area of ocean, the variation in specific gravity 
values is not large; it is near the average specific gravity of quartz and 
the carbonate type minerals. The determination of the unidimensional 
compressive strength, q,, besides the previously discussed soil mechanics 
regulation for carrying out testing, is according to the Ministry of Geology 
and Mineral Resources, System Unit 9 Rock and Acoustics Group's methods 
regarding the handling of test data. Table 2 includes important parameters 
generally used to establish the ocean floor sediment elastic model which is 
usually used. Table 3 lists the sound velocity values of several of the 
present experimental ocean area's representative test sites as determined by 
reflectance methods2). Sound velocity data from the other test sites can be 
found in Table 2 of Reference [11]. 





2) In this experimental ocean area, sound velocities below the ocean floor 
surface were all determined from reflectance data. 
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Table 3. Some Representative Test Site Sound Velocity Data From the Experi- 
mental Ocean Area (based on Reference [11]) 
Test Site 
Designation T S F V 
Water Depth (m) 20 39.5 60 69 
Month of 
Observation Oct 1980 Oct 1980 Sep 1981 Mar 1980 
Sound Velocity in 
Bottom Layer of 
Ocean Water, co(m/s) 1539.7 1540.3 1538.7 1530 
Sound Velocity in Depth; Sound | Depth] Sound | Depth | Sound Depth] Sound 
Layers of Ocean (m) speed | (m) speed| (m) | speed| (m) | speed 
Floor Sediment, (m/s) (m/s) (m/s) (m/s) 
cj (m/s) 
(depth is figured 9 1514 7 1520 ' 1530 1520 
from the ocean 1.94 ~1685 9-2 ~1575| 2.35 1.3 1543 
RAOGE SaEtere) 98 30 3.4 | 1663] 1's] 1551 
>98 | 4800 | >30] 1705] >3.4 | 1663] 2/3] 1569 
2\7 ~ 1585 
' 
36 
>36 1701 





























IIl. 


Mechanical Factor, q,,, Which Influences the Velocity of Sound 


R.S. Anderson!5] has extensively investigated each physical-chemical property 


of sediment which influences sound velocity. 


These include the specific 


gravity, G, the carbon content, C, the porosity, n, the water content, w, 
the porosity ratio, e, the average particle diameter, $, the particle group 
kurtosis, K, the skewness, SK, the weight of moisture, etc., but, except for 
porosity, none has given very satisfying results. 
data calculated from empirical formulae dealing with the porosity-sound 
velocity relationship (for example, the double parameter equations, f(c,n) = 0, 
of Anderson and Hamilton, etc.), still there is a relatively large difference 
in values with the actual sound velocity measurements due to non-linear 


relationships. 
et al. in 1983(11 


Merely according to the 


This was already pointed out by Zhou Zhiyu [0719 1807 1946] 
In addition, the results we obtain in the laboratory by 


repeatedly applying acoustical rock parameter determining instrumentation to 
samples, clearly show that the difference in sound velocity with the same 
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sample (horizontal and vertical directions) at its maximum could reach 150 to 
200 meters/second; this cannot be explained using porosity alone3). Thus, 

it is clear that, besides porosity, trere are still more independent parame- 
ters which influence sound 


~~ +e 4 & F © 


Based on previous analysis, we differed from Anderson and adopted other 
approaches to a solution to the problem. We took the soil mechanical 
parameters of sediment, including the flow limit, Wr, the plasticity limic, 
Wp» the plasticity index, Ip, and the unidimensional compressive strength, 
du» etc., along with the porosity parameter, n, and carried out regression 
analysis on the relationships to sound velocity, c, via a computer. Results 
clearly show that taking the sediment compressive strength, qy, as a supple- 
mentary factor when using the porosity, n, to calculate sound velocity, 
compared to merely using the porosity factor, there is a better degree of 
correlation. For example, in comparison with Anderson's empirical formula! 
which is based on just ocean sediment sound velocity and porosity, when ay is 
added as a third parameter and applied to calculation of the sound velocity 
in this ocean area, the standard deviation is lowered to 14.9 m/s from the 

35 m/s of the Anderson formula; relative error is less than 9.95 percent. 
Then the three parameter regression formulas we obtain from our work in this 
ocean area has the following form: 


/? 


f(c,m,g.) =C— 4195490 5833m—0 7 695n° +94, 908g, — 64, 460395 = (4) 


Our calculation model is shown in Figure 4. Taking experimentai ocean area 
test site V as an example, a comparison of calculations of sediment sound 
velocity using Equation (4) and other often used formulae can be inspected 
in Table 4. 


The implications of Equation (4) are such that in practice there is the 
possibility for the statistical characteristics of microscopic pore struc- 
ture to be reflected in the macroscopic whole, and supplementary revised 
results of ocean floor sediment sound velocity calculations will be 
obtained. qy is just this kind of mechanical factor which can reflect 
microscopic pore structure as well as the nature of the granular buildup 
and the comprehensive effects of the state of mechanical coupling, when 
used as a supplementary term in calculations of the sound velocity, c. 


In order to further examine the utility of Equation (4), we used it to 
calculate data from layers from several other positions in this ocean area. 





3) Liang Yuanbo and Lu Bo, “Acoustical and Mechanical Properties of Shallow 
Ocean Floor Sediment,” 1983. 
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Test Site V Media Parameters 




















Average 
Layer particle 
thickness Porosity diameter 
Layer # (cm Sediment type (%) (mm) 
I 34 sticky earth 57.0 0.0810 
II 36 powdery sand 56.6 0.0545 
III 24 sticky earth 59.3 0.0560 
VI 36 sticky earth 48 .6 0.0680 
Ocean water 
Sound velocity of bottom layer 
Ocean floor of ocean water, cg = 1530 a/s 
0 Area 1 
(cm) sedimentary sound velocity, cj = 1520 a/s 
I average porosity, n = 55.4% 
Wh Ke shenewesetse ews unidimensional compressive strength, 
II qy * 0.078 kg/cm2 
ME RURESHR ASS S SEES calculation formula: cc. = 4195 - 90.5833n 
Il + 0.7695n2 - 94.6968q,, + 64.4603q,2 
OF cccccccccccccccce 
VI 
130 Area 2 





| cp = 1543 o/s 


Figure 4. Media Parameter Model 
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Table 4. Comparison of Various Caiculation Results on Sedimentary Layers in 
the Experimental Ocean Area 
On-site Calcu- 
Experi- Unidimen- lated 
mental Sional Results, 
Sound Compressive Sound 
Test Velocity Porosicy Strength Velocity Deviation* 
Site c)(m/s) n (%) q,, (kg/cm-) Formula cc (m/s) \c (m/s) 
l. Anderson 
Formulal®l; 
c=2506-27 .58n 1551.4 +31.% 
+0 .1868n¢ 
2. Hamilton 
Formulal8); 
c=2455.9-21.716n 1639.5 +119.5 
+0.126n2 
\ 1520 55.4 0.078 
3. Formula cited 
in Reference [2]: 
c#=2475.5-21.764n 1647.3 +127.3 
+0.123n¢ 
4. Formula (4) of 
this work: 
c=4195-90.5833n 1531.3 +11.3 
+0.7695n2- 
94.6968q, 
+64 .4603q,,2 
the = C, - Cc] : 
Results were as follows (same units as in Table 4): 
Test Site Layer (cn) n qu _ Cc. C] \c 
F 43 9.1 0.151 1516 1530 14 
L 60 62.2 0.451 1507 1520 13 
U 40 62.4 0.754 1504 1514 10 


within the range of the standard 


c is 








leviation. 








IV. Some Discussion 


The physical significance of Equation (4) is easy to understand because taking 
qu as the upper limit of a sediment's resistance to vertical pressure, it is 
certainly one of the intrinsic physical-mechanical properties (or charac- 
teristics) of natural sediments; it is, due to the granular composition of 
the porous media, a determining factor of the pattern of granular buildup 
and the state of mechanicai coupling. Then, speaking of relatively broad 
diameter particles on powdery sand, the important effective elements are the 
internal friction of the solid particles, and the cementing effect. As for 
the pattern of granular buildup, aithough Graton and Fraser|12] analyzed 
buildup in simple spherical systems and found that there are six common and 
important patterns, and that sediment is influenced by density and unit 
weight, natural sediments rarely have isometric, homogeneous, spherical 
buildup, leading to a situation in which there is no way to obtain a full 
set of factors through concrete testing (see plates 2, 3 and 4), thus, they 
cannot simply be accounted for using density or unit weight. Additionally, 
due to the cementing activity in back of the sediment (plates 5-8), one is 
also unable to use suspended solids equations for mineral particles in sea 
water to describe sedimentation (Reference [16], pp 7-11). 


We have taken another look at what happens in the situation where we merely 
use q, and don't consider the porosity. Gen An [2704 1489] and others have 
dealt with rock cores from Hokaido, Japan as well as rock experiments carried 
out in several of our large test areas by China's Ministry of Geology and 
Mineral Resources, System 9 Unit Rocks and Acoustics Group. Their results 
are as follows: 


lg(q.) =2.671lg(c,} +! .09 
Or, . go =12,3c;°°" (5) 


which omits the effect of the rock porosity on sound velocity. This type of 
double parameter formula (f(c,q,,) = 0) deals with this problem from a differ- 
ent approach, but its results are applicable only in cases of the lowest 
porosity (for example, n < 10%), since in this case the sound waves are 
basically propagating through a solid framework. Morecver, the few empirical 
formulae discussed above apply to test results of igneous rocks and meta- 
morphic rock layers. But, speaking in terms of the water-saturated porous 
media of contemporary ocean sedimentation, the porosity, n, cannot be 
neglected. This point has already been amply expounded upon by Biot 6], 


In common with any of the empirical formulae, Equation (4) was established 
from test data from specific experimental ocean areas and its coefficients 

in different ocean areas can have some variation. Our work, and the above 
analysis, just point out une way to deal with this; namely qy is an effective 
factor which, together with the porosity, n, conditions the velocity of sound 
propagation in ocean waters. It naturally includes the macroscopic effect 

of the pore structure of the porous media. The results of our regression 
analysis of this area's sediment test data clearly show that this is the 


case. 
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VY. Conclusion 


1. Adopting the mechanical factor, q,,, comprehensively expresses the overall 
effects of granular buildup of ocean floor sediment and pore structure and 
can serve as an additional corrective factor in calculations of sedimentary 
layer sound velocity. It raises the calculation precision of the often used 
formulae which don't consider pore structure, and at the same time, it 
eliminates the necessity of carrying out statistical measurements of micro- 
scopic pores for correcting sound velocity calculations. 


2. This work and related analyses make clear that the use of ocean floor 
acoustical characteristics to distinguish the engineering and mechanical 
properties of sedimentary layers is not only possible, but also there are 
broad areas which need investigation and research; for example, investigation 
of the relationships among the shearing strength, S, the Adebao [7093 1779 
1027] liquid limit, W,, acoustic impedance, Z, etc. As to the problem of 
sedimentary shear wave propagation and shear modulus, yu, we will discuss that 
further in a future paper. In essence, these are established through 
acoustical methods together with sediment mechanical properties and work 
dealing with relationships among geological factors. With this aspect of 
investigative research we can hope to more perfectly develop the technology 
for remote sound sensing and telemetry of ocean floor material. 


The on-site sound velocity information used in this report was provided by 
the lst Laboratory of the Ocean Floor Group (Comrade Zhou Zhiyu [0719 1807 
1946] and others), a cooperative unit of the Institute of Acoustics, Chinese 
Academy of Sciences. Also work on the electron microscope observations and 
picture taking was received through the enthusiastic cooperation of the 
Electron Microscope Laboratory, Zhongshan University; for this, we express 
our deep appreciation. 


PLATE CAPTIONS 
Electron micrographs of ocean floor sediment structure: 
1. Structure of near spherical and equal granular buildup (consisting of 
pyrite particles. South China Sea, northwest section continental shelf, 


ocean floor area #551, magnified 30,000 times. 


2. Structure of granular buildup of natural sediment. [bid., /262, magni- 
fied 5,000 times. 


3. Structure of sedimentary pores having definite orientation. Ibid., #336, 
magnified 30,000 times. 


4. Microscopic structure of plate-like bioclastic sediment. Ibid., #328, 
magnified 7,500 times. 


5. Cementation of sediment particles (cementing of aragonite needles). Ibid., 
#686, magnified 8,000 times. 
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Microscopic structure of aragonite needle cementite. Ibid., #683, magni- 
fied 7,500 times. 


Skeletal structure of network shape formed by aragonite needle cementite. 
Ibid., #352, magnified 30,000 times. 


Filling and proliferating conditions of aragonite needle cementite among 
granular pores. Ibid., #361, magnified 15,000 times. 
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PHYSICAL SCIENCES 


COMPARISON OF DEEPWATER WIND-WAVE PREDICTION MODELS 


Qingdao SHANDONG HAIYANG XUEYUAN XUEBAO [JOURNAL OF SHANDONG COLLEGE OF 
OCEANOLOGY] in Chinese Vol 15 No 3, 15 Sep 85 pp 1-11 


[Article by Wang Bingxiang [3076 3521 4382], Zhang Jiuzi [1728 1432 5417], 
Department of Physical Oceanology and Marine Meteorology, and Wu Guiqiu 
[2976 2710 4428], First Institute of Oceanography, National Bureau of 
Oceanography: "Comparison of Deepwater Wind-Wave Prediction Models"; 
paper received 5 May 1984] 


[Text] English Abstract: Presently, there are about 36 models used in deep- 
water wind-wave forecasting by incomplete statistics. To select from among 
these models 12 that are more important and currently employed are con- 
firmed by field observations, and are compared with each other. The calcula- 
tions based on the 12 have been compared with wave data proffered by Iu.M. 
Krylov (1958) and Bao Qiangsheng et al. (1983); so far as the cases 

examined are concerned, the accuracy of the results given by Wang's model 
are superior to those of the other 11 models. In order to investigate 
typhoon waves at station "Xiaomaidao" at 1800 hours 30 August 1939, the 
above approaches to wave prediction are used to hindcast again; the wave 
height on a basis of Wang's model is found to be more consistent with a 
50-year return period wave at the wave station than those found for the 
existing models. According to the prediction technique of hurricane waves 
the wave height and period of the most probable first highest wave within 
the area of the typhoon are estimated at about 39 meters and 24 seconds 
respectively. 


Ocean wave wave height and period are two important parameters in marine 
engineering site selection and structure design. There are generally two 
applications of these two wave factors obtained for a designated region: 

One is if there are already wave observation materials, then spectrum or non- 
spectrum methods can be relied on; the other is if there are no materials, 
then with the help of wave materials for neighboring marine regions, they 
can be deduced using refraction-diffraction theory or, after determining 

the wind factors from historical weather maps, they can be hindcast using 
wave forecast methods. We followed the second method in the second cate- 
gory to estimate the marine typhoon wave factors. 
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Since Sverdrup and Munk first published their method for forecasting effec- 
tive wave in 19741, many methods have been proposed. Counting only those in 
references 2-10 there are 36. Although there are so many forecasting models, 
they can be reduced to two categories: quasi-experimental and purely 
experimental. We selected 12 of the more important ones which are frequently 
used, for verification and comparison and carried out a wind-wave hindcast 
based on a strong typhoon in the history of the Huanghai. 


For the verification which was carried out on each method using the materials 
published in Table 2 of reference 6 and Table 1 of reference 11, the method 
of reference 4 provided better results than the other methods. 


The methods were applied to hindcasting typhoon waves at Xiaomaidao at 1800 
hours on 30 August 1939; based on the results of the method in reference 4, 
conformity to the observation values was better than the other methods. 

As concerns the results obtained using the hurricane wave forecast methods, 
the maximum probable values of the wave height and wave period within the 
typhoon area are 39.3 meters and 24.19 seconds, respectively. 


I. Verification and Comparison of Models 


1. The models were verified on the basis of the data provided in Table 2 of 
reference 6 and Table 1 of reference ll. 


Apart from the eight models in references 3, 4, and 6, 4 of the 12 models 
used were newly selected for uniform units, and after conversion were 


(1) G.B. Rzheplinskii (1972)!?: 
(Hy )e= 1,.7872K 10° etsy! 


Ts = 0,1126+U%7% Sey! 75 
(Hey): = 6.6911 10 4eppi 97 lenges - 


Tr = 6.1437 10°20 U?“4eoqi 4 , 
(2) B.W.Wilson, WV (1965)''°:; 
(Hex )a= 3.7277 10 2*U2eC1—(1 + 1,2528K 10 Ze fox! 2) Y 
(3) 


Te = 0.7590°UC1—( 1 1.7126 10 -*U 2 Sex! /*) 8 
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(3) D.J.T.Carter(1982):''°, 


J Hews 6.2821 x 10° 4+yex!“- 


,tes— —_ } (4) 
Tx = 5.5267 1072 +Yete yo. * 





—Yes@—¢ 
\ J Hawda= 3.0225 x 10° ?+y- 
“° “T= 0.566 } (5) 
D : U 
(Hex): =5,1292x 107569’ Tog577 
7 res on (6) 
4 T,. =1,2534~x 10°20 [47 749377 j 
— No a 
\ J Haw) = 2.9250x 10°*-y2.. 
= , (7) 
Te = 0,566°U é 


in which 


D 1>33.3713-y %seg%e7 
@ x<2,32x 10?-y? 
@ t<7.2306x 10?+y 


(4) The method in reference 2 (abbreviated QQQ) (1980): 


(Hoy )e= 3.4786 10 7*U?+C1— exp( — 2.5057 10 2+ tex! 2) ) 
_ . (8) 
Te = 3.645 Hox 


The units U (wind velocity), x (wind area), and t (wind time) in equations 
(1)-(8) are, respectively, m/sec-l, meters, and seconds; the Hsy (accumula- 
tion rate is 5 percent of wave height) and T (average period) units are 
meters and seconds, respectively. 


The wave data in Table 1 of reference 11 was used in references 3, 4, 6, and 
13; the wind factors ,iven in this table were substituted in the relevant 
forms in equations (1)-(8). The results of the wave factors calculated are 
given in Table l. 


From the figures in Tables 2 and 3 we know that the wave height given by 
the method in reference 4 is better than the other methods. 


Here we should explain that conversion between the wave heights was induced 
according to theoretical resultsl5, i.e., 


H,., - 1.95H 
HW, = 2,032 ) 
To (9) 
Hi,3= 1,598H 
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1. Depth (meters) 5. Wave height (H5z7) meters 
2. Wind velocity (m/s) 6. Period, seconds 

3. Wind-time (hours) 7. Observations 

4. Wind area (km) 


Note to Table 1: The letters SM, B], K, PNJ, WEN, WANG, B2, NOR, P, C, WIL 
and QQQ in Tables 1, 2, and 3 represent the methods of Sverdrup-Munk (1947), 
Bretschneider (1952), Krylov (1958), Pierson-Neumann-Japnes (1955), Wen 
Shengchang (1960), Wang Bingxiang (1964), Bretschneider (1970), norm (1978), 
Rzheplinskii (1972), Carter (1982), Wilson (1965), and the Shandong College 
of Oceanography-First Institute of Oceanography-Institute of Oceanography 

of the Academy of Sciences (1980), respectively. When these letters appear 
later in this paper, they carry the same meaning. 


Table 2. Maximum Error of Wave Height Calculated According to Table 1* 











; cae VgmneeRe 
= | NO 
SM 108.6 | i 
B: 82.9 | 41 
K 73.6 4! 
PNJ 75.0 | 18 
WEN 41.1 | ? 
WANG 35.0 41 
Bs 61,0 22 
NOR $1.0 | 40 
P 61.4 i 
Cc 50.0 | 
WIL 133.2 40 
QQQ 146.3 | 40 











*SM, B,, K, PNJ, WEN, WANG and B?, 
NOR are taken from references 4 and 6. 


Key: 
Ll. Method 
Maximum wave height error 
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Table 3. Maximum Error of Wave Height of the Various Methods Calculated in 
Table 2 According to Reference 6% 














, ss ‘-) ganache 
* | NO 
SM 100.0 | 12 
B; 98° | “ 
K $2.0 | t 
PNJ 105.0 : 
WEN 46.0 | 3 
WANG 32.0 5 
Bs $$.0 2 
NOR $1.0 12 
P 61.8 : 
c 42.3 : 
WIL 81.7 1 
QQQ 93.2 | : 











*Cited according to reference 6 
except for P, C, WIL, and QQQ. 


1. Method 
2. Maximum wave height error 


The relationship between periods is according to experimental equation! 


T, = 1,156 T (10) 
3 


This value is almost equivalent to the 1.15 of reference 17. 


2. It was difficult to carry out a systematic comparison of the models. 

The main difficulty was that of the 12 methods to be compared in the paper, 
if we relied on the wind-wave dimensions for their growth stage we were 
restricted to three essential divisions of wind, the three were: 1, siml- 
taneously relative wind area and wind time growth4; 2, only relative to wind 
area growth2»10; 3, growth relative to either the wind area or wind time. 

To overcome this difficulty and get on with the comparison, first of all, 

we handled it as reference 3 in comparing this item, and selected only one 
representative wind velocity, for example, 15 m/sec-l; next, since there 
were two in the paper which were only relative to wind area growth, and 
therefore for wave theory of growth relative to their wind area or wind 
time, we took only the expression for the former; for processing reference 4, 
two wind times were provided: one was the wind time (3.42 x 10? sec) 
necessary for the wind wave to achieve full growth at the corresponding 
wave velocity of 15 m/sec~!, and the other was 102 sec, which was changed 


into only relative to wind area growth. 
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In Figures 1 and 2, the series marked WEN, SM, By, K, and PNJ reproduce 


Figures /.7-3 in reference 15, the rest are calculated and drawn using 
various methods. 
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Figure 1. Comparison of Wave Height of Various Methods With Growth of Wind 


Figure 2. 
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Comparison of Wave Period of the Various Methods With Growth of 


Wind Area 


From the curve in Figure 1 the following characteristics can be seen: 


K curve from beginning to end is within the upper and lower limits of the 
Wang curve; 2) the Carter and PNJ curves are both below the lower limits of 


the Wang curve when the wind region length is smaller than 10° cm and 
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- x 10° cm, respectively; 3) the SM curve first goes 30 km beyond the upper 
limit of the Wang curve, and continues to go beyond it, next are the Bj, 
Wilson, Rzheplinskii, NOR, B?, Carter, and PNJ curves, and finally 300 km 


1, se : so «Ft Lyra” ‘ 
Jc} ond 1C 1s che wl. curve. 


Because the Lines in the figure are so concentrated, there isn't any way to 
draw the b> curve, therefore it is explained below: 1) within the 3-200 km 
range, Bz is slightly higher than WEN and slightly lower than K; 2) from 
90° km it almost overlaps WEN and is higher than K; 3) in between 3 and 
2°00 au, is lower than B}. 


The curves in Figure 2 reveal: 1) By initially goes beyond the upper limit 
of the Wang curve, those which continue to go beyond it next are NOR, SM, 
zheplinskii and Wilson; 2) in the wind region under discussion, WEN and K 
hroughout are located between the upper and lower limits of the Wang curve; 
3) SM and Wilson are below the lower limit of the Wang curve when the wind 
area is smaller than 6 km and 10 km respectively. 


,m we 


Three interesting relationships must be explained here: 1) At 4x103<x<104 cm, 
the Carter and Wilson period curves are both almost congruent; at 2x10%<x<10/ 
cm, the former is lower than the latter; at 107<x<108 cm, it is the reverse; 
when X>108 cm, the former is again lower than the latter. 2) At X<440 km, 


NOR is higher than SM and Rzheplinskii; at X>400 km, it is the opposite. 
3) At 3<X<60 km, B2>SM; at X>60 km, K<B2<SM; By is always lower than B). 


From Figure 1 and Figure 2 it can be clearly seen that wind area and wind 
time have an influence on the wave height and period computed on the basis 
of the method in reference 4; when the wind area is short, wind-wave dimen- 
sions increase very rapidly with the wind time, and when it is long, the 
rule played by the wind time gradually diminishes until X = 3240 km, the 
wind time effect becomes so weak that even the wind height and period at 

th = 95 hours and t = 102 seconds, are respectively 2.68 m and 2.00 m and 
8.22 seconds and 7.10 seconds. 


Il. Hindcast of Huang Hai Typhoon Waves 


In summer and autumn China's Huang Hai is occasionally hit by typhoons giving 
rise to typhoon waves, which have a catastrophic impact on China's coast!8, 


survey of historical materials on typhoons reveals that in the 84 years 
of this century the one whose impact on the southern coast of Shandong was 
most severe was the typhoon that went through the Huang Hai on 28-31 August 
1939, It goes without saying that the dimensions of the typhoon waves which 
accompanied it also should be the largest. So, how high and how long were 
its wave height and period? On the one hand, this paper used the 12 models 
verified and compared in the previous section to hindcast the typhoon wav 
at Xiaomaidao at 1800 hours on 30 August 1939, and on the other hand rel 

m wave factors provided by the hurricane wave forecast methods!9 for the 
first largest wave likely within the typhoon area, we attempted to give a 
probable answer to this question. 











The present problem is how to determine the wind area, wind velocity, and 
wind time required by wind-wave forecasts. 


Since the wind velocity materials of the surface temperature map of 1800 
hours on 30 August 1939 given in Figure 3 is far from able to Satisfy the 
demands of wind-wave forecasting, we followed the read-outs of the wind 
velocity curve of the Qingdao meteorological station and after calculation 


treated the values obtained as the wind velocity required. 











Figure 3. Weather Map of 1800 Hours Figure 4. Chart of Meteorological 
(Beijing time) 30 August Changes at Qingdao 
1939 30 August-1 September 


1939 (cited from Figure 3 
in reference 18) 


From the curve in Figure 4 we read the average wind velocity at 1200 hours 
and 1800 hours on the 30th as U; = 17.5 m/sec-! and U2 = 27 m/sec7!, then 
uSing the formula for when the wind velocity continues to increase we solved 
for the representative value between 1200 and 180029: 


U(Z})] = 0.2U; + 0.8U2 = 25.1 m/sec™! (11) 
which is 0.6 m/sec! greater than the lower limit for a force 10 wind (24.5 
m/sec~!). Experience tells us that the wind velocity over the sea is 
greater by one or two grades than over the land. Thus, we took its upper 
limit, i.e., 2 grades, because from Figure 3 the wind area under considera- 
tion was situated on the right side of the typhoon's direction of aerneees 
which was the area where the wind velocity was greater than its left side 9, 
Therefore, the wind at sea corresponding to the Qingdao Guanxiangshang wind 
velocity meter, a point 103.8 meters above sea level (Z}) the wind was the 
lower limit (32.7 m/sec~!) of force 12 plus 0.6 m/sec@!, i.e., U(Z)) = 33.3 


m/sec=! ° 


~~ 


According to the vertical distribution of wind velocity, we found the wind 


velocity at a point (Z?) 10 meters above the surface of the seal8; 








UCZ:) = U(Zi )-| EE Hees = 26,54 n/sec~1 (12) 


in which Zp = 0.001 is the Ellison's observation value at the sea surface. 


The wind area relies on the idea of the most favorable line wind areal’, and 
is determined by reference to the wave propagation direction, i.e., in 
Figure 3, AB = 350 nautical miles. 


Determination of the wind time will involve equivalent wind time te, for the 
same wave height and wave velocity, different methods will give individual 
values, and thus are related to the wind-wave dimensions sought. For pur- 
poses of simplification, this paper takes the norm as the standard. Assuming 
that at 1200 hours on 30 August the typhoon had not yet reached the area 
under discussion but already had waves with wave height of 5 m8. This value 
together with wind velocity U(Z2) can be found from the norm te = 2.2 hours, 
thus we obtain the forecast of 30 August's wind-wave wave time as 8.2 hours. 


The wind factors obtained were substituted in the 12 models in the previous 
section to hindcast the wind-waves at Xiaomaidac at 1800 hours on 30 August 
1939. The results are given in Table 4. 


Table 4. Comparison of Various Methods for Wind-Wave* 
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*These data are cited from reference 21 Xiaomaidao 50-year return 
period design wave factors. 


Key: 
1. Time: August 1939 5. Wave height (H)/;9) meters 
2. Wind velocity: (m/sec~!) 6. Period (T) seconds 
3. Wind-time: (hours) 7. Percentage error (7%) 
4. Wind area (km) 


If we make the Xiaomaidao 50-year return period design wave factor as the 
criterion, then from the figures in Table 4 we can tell that: 1) In terms 
of percent error, the errors of the various methods are smaller than the 

corresponding methods in Table 2 and Table 3, this seems to mean that the 
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method of determining wind factors is rational and thus the values obtained 
conform to reality. 2) The wave height values provided using the model in 
reference 4 are better than the other models, but the best period values are 
at the same time provided by the norm20 and Wilson methodl0. 3) The period 
error is smaller than the wave heights. 


The model for forecasting deepwater hurricane wavel?: 








R 
P,—P, = “— 
P.—-P, ° (13) 
User = 0,868(73*AP®*°— R*(0,575f)) (14) 
Ur = 0,865U esx + 0.5VF (15) 
B-AP 
H1=16.5¢ 1 -[1+—2-208aVe] 
/ Ue (16) 
T 1 =2,13/ Hy 
Tv 7 (17) 


in which: f (Coriolis parameter) = 2w:sind; if ¢ is the latitude of the 
center of the typhoon, 33.4°N, and w (the earth's angular velocity)= 21/24 
is substituted, we get f = -0.285 arc*hours-l 


a is a parameter which is determined by the forward velocity of a hurricane 
and the growth of the effective wind area, and for a hurricane moving slowly, 
it can be assumed to be l. 


constant pressure py, = 29.92 hours!9, 
po and p are both pressure, units both are hours. 
Ap = Pn-PO- 
units of effective wave height (185) and period (15) are feet and 
seconds. 


If we substitute the data read from the surface weather map at 1800 hours on 
30 August 1939: 


the typhoon center's distance r from Qingdao = 360 km; 


the atmospheric pressure P¢ at a point r distant from the typhoon center, 1i.e., 
the Qingdao Station, = 995.9mb* and the data cited in reference 18: 


the atmospheric pressure pg in the center of the typhoon = 959.9mb; 
the center of the typhoor advanced at the average velocity 
VF(=8.095 knots) 





*Document [18] misread the 30 August 1500 hours weather map time as 1800 
hours, thus the corresponding errors are 434 km and 998.6mb. 
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in equations (13)-(17), then we get: 


maximum wind velocity radius R = 76.31 nautical miles; 

maximum wind velocity gradient Up;, at a point 30 feet above the sea 
surface = 68.43 knots; 

for the moving wind area, the maximum continuous wind velocity Up at 
a point 30 feet above the sea surface = 63.24 knots; 

the effective wave height and period at the place of maximum wind 
velocity are, respectively, 


H 20.19 meters (18) 


Wo | 


La 
Wo] 


17.34 seconds (19) 


Assuming that according to reference 19 the most probable largest sea wave 
is determined by the number N of waves; and N is determined by the length 
of the hurricane's cross section in a zhundingchang [0402 1353 1603] state 
and its forward velocity. It has already been shown that the largest sea 
waves appeared in a range equivalent in length to the radius of the maximum 
velocity, and the time needed for this length to go through a typical point: 


= 33936 sec (20) 


N= ——=1957 


Ti+3 (21) 





the probable maximum wave height, is calculated according to 


— 

H, = 0.707 x Hifin(-®) (22) 
b) n 

take n = 1, then we get the probable first maximum wave height 


H; = 39.30 meters (23) 


lf we substitute this value in equation (17), then we get the corresponding 


T; = 24.19 seconds (24) 


Conclusion and Discussion 


According to incomplete figures, up to the present there have been 36 deep- 
water wind-wave forecasting models proposed. We have selected the more 
important and often used 12 for verification and comparison and carried out 
wind-wave hindcasting for a powerful typhoon that hit Huang Hai, and obtained 
the following: 
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(1) In terms of testing the models using the observation data provided by 
references 6 and 11, the wave height value calculated according to the model 
in reference 4 was better than that of the other models, the wave height 
percentage errors were 35 percent and 32 percent, respectively. 


(2) We carried out a hindcast of waves in a typhoon in Huang Hai, and using 
the design wave factors of a 50-year return period wave at Xiaomaidao as the 
norm, then the method in reference 4 gave better wave height values than the 
other methods, and the best period value was provided by the norms and 
Wilson IV method simultaneously. 


(3) As for which ratio should be selected for T1/3/T is of general signifi- 
cance and will have to be explored further in the future. 


(4) The values obtained by converting the materials read from the wind 
velocity record curves of a single station on the seacoast to serve as wind 
velocity in a wind area have more significance and also can only be deter- 
mined after systematic exploration in the future. 


(5) As for the examples examined, although among the 12 models selected the 
method based on reference 4 gave better results, before reaching any final 
conclusions, more extensive and systematic verification will still be 
required. 
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PHYSICAL SCIENCES 


METHODS OF FREQUENCY SPECTRUM ESTIMATION OF SEA WATER COMPARED 


Qingdao SHANDONG HAIYANG XUEYUAN XUEBAO [JOURNAL OF SHANDONG COLLEGE OF 
OCEANOLOGY] in Chinese Vol 15 No 3, 15 Sep 85 pp 12-21 


[Article by Jiang Decai [5592 1795 2088], Zhang Dacuo [1728 1129 6934], 
Department of Physical Oceanology and Marine Meteorology, and Gao Xinsheng 
(7559 2450 3932], Computer Lab: "The Comparison of the Maximum Entropy 
Method With Conventional One for the Frequency Spectrum Estimation of Sea 
Wave"; paper received 31 May 1984} 


[Text] English Abstract: The principle of the MEM (Maximum Entropy Method) 
is briefly introduced in this paper. The determination of the length of 
prediction filter is discussed in the case of the MEM being used for esti- 
mating the spectrum of sea wave (wind wave, mixed wave) profile. The 
quality of the spectrum estimation with MEM and conventional methods is 

also dealt with under the condition of the same sampling interval and 

the same sample capacity, it is shown that the MEM has higher resolution, 
especially for the short series of wave profile records, it is much better 
than conventional methods. 


1. Statement of the Problem 


The sea wave power spectrum is a basic characteristic parameter in describ- 
ing the random processes of sea wave wave surface, it can study the distribu- 
tion of sea wave energy relative to frequency, as well as study the statistical 
laws manifested by sea wave wave surface. Thus, it is very important for 
accurately determining the sea wave frequency spectrum. The correlation and 
period map methods are the common methods currently used for estimating the 
sea wave frequency spectrum. They both introduce the subjective will of some 
hypotheses and certain weight functions. The Maximum Entropy Method (MEM) is 
a method in which the materials adapt to spectrum estimation. Its window 
functions are not fixed, it adapts to one frequency and then another, the 
estimated spectrum values at one frequency are rarely disturbed by variances 
of other frequencies!. In 1972, Ulrych2 applied Burg? MEM to estimate 
frequency spectra, demonstrating the use of short record materials in 
spectrum estimation, it had high resolution, and was superior to the common 
correlation and period map methods. This paper applies the MEM to processing 
random data of sea wave surface and compares it with the traditional methods 
of estimating sea wave frequency spectra in the search for a aore effective 
method of estimating sea wave spectra. 
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[The MEM has been applied to many geophysical areas with great results. 


There have not been many application to sea wave frequency spectrum estima- 
tion. Takahashi Tomoharu et al.* (1977) used MEM to estimate the frequenc\ 
spectrum of nonstationary ocean wave surfaces and obtained preliminary 


results; Tsurutani Ryuichi et al.? (1977) used MEM for frequency estimating 
of wave surface records of wind-wave trough measurements; Sverre et al.® 
used MEM to estimate ocean wave surface frequency spectra and discussed the 
relationship between the spectrum parameters (spectrum moment) and the final 
prediction error filter length. The work of the above-mentioned scholars 
showed that using MEM to estimate sea wave frequencies had the advantage of 
high resolution, but there is still one common problem, i.e., determining 
the length of the prediction error filter in MEM requires combining specific 
circumstances of the wave surface (wind-wave, swell, mixed waves). This is 
because if the filter length is too short it gives a smoothed spectrum, and 
if it is too long, it also gives secondary peaks formed by false signals. 
This paper will combine some wave surface records measured by different 
instruments in China and use the MEM and traditional methods to estimate 
spectra, then compare the spectral forms, spectral moment, and wave factors; 
discuss the method for determining the final predictioa filter length--the 


} 


Akaike fixed stage method and capacity ratio fixed stage method to derive a 


more effective method of estimating sea wave frequen pectra. 


il. Basic Principles of MEM 

MEM estimates directly from the data a minimum phase prediction error 

filter, and solves the input power spectrum by using the output power 

spectrum P(f) of the prediction error filter and the prediction error 
° ’ 

filter spectrum |H(f)/¢ 


S(f)=, P(f) 

Hit f) (1) 
MEM is a method of extrapolating the auto-correlation function to the outside 
of the maximum time lag, its basic idea is to find the entropy maximization 
at the time of each extrapolated auto-correlation function. It is not 
restricted by the length itself of stationary random process information, it 
relies on the predictability of the sampling of both ends for its better 
resolution, that is, when the entropy is maximized, the auto-correlation is 
extrapolated. 


The relationship of entropy and the stationary Gaussian process spectral 
density derived according to Smyliel! (1973) 


rh 


(2) 


atk 
if P| _ Mee tees const 


in which H is the entropy of the stationary Gaussian process and ty Is 


Nyquist frequency. 











According to the Wiener-Khintchine relationship 


iy 
| _ S(OZMd t= R(kAt) (3) 


in which Z = exp(i2nfAt), At is the uniform sampling step length, R(kAt) is 
the auto-correlation function which has a time lag (kit). The MEM ends up 
as a variation problem of solving equation (2) relying on equation (3) as 
the limiting condition, and the solution to the variation problem is 


s(f=(Z 4,2") (4) 


in which Ay is the Lagrange multiplier, determined by equation (3). The 
ocean power spectrum (frequency spectrum) is a real number and non-negative, 
thus equation (4) can be rewritten 


Sire -——FelSte. ___ (5) 


i= ¥ ae2Z" 
i! 


in which (i- E a2") is the minimum phase, Pp, is the average output power 
kel 


of the m+ 1 prediction error filter. Rewriting equation (3) as 


S(P=5 —— x R(kAf)Z* (6) 
a e-e 


comparing the equal coefficients of Z in equations (5) and (6), and taking 
into account that R(kAt) is an even function, we get 





(ak l { » | 
| pidsaliiceatie aii =| ' (7) 


the first matrix on the left is the Toeplitz matrix. The minimum phase 
filter coefficient am, and the average output power P_, are determined by 
equation (7). 


Suppose that the data of an actual sea wave surface event is already known 
and is expressed as €], C25 +++» Gn» the forward and backward prediction 
error filters designed by Burg’ (1967) are the coefficients an k determined 
by the prediction error square average (output power) minimization. The 
output power P) of the two prediction error filters (1,a),)) wave surface 
sequence (;(i = 1,2,...,n) through the forward and backward prediction error 
filters 
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‘ Staay = (¢ te Gist + (Gio — a O27 (8) 
-1 


in which (C;-a) »ASi+1) is the ¢; backward prediction error, (044)-a1,1%4) 
is the ¢; forwar prediction error, and if we wish to minimize the average 
output power, then 2P)/2a),1 = 0, we have 


ry s | 
aint = 2H Cini / FCCP +E, D (9) 


From equation (9) we can see that a)_/}<l is the minimum phase prediction 
error filter. 


From equation (7) we get 


Ri = a151 Re 
Pi = (l—aj.1) Ro (10) 
R= ¥ei 
- s 1 
and applying the Levinson recursion formla 
Bert =e ict ~ Aasade isa-t (11) 


which means that the filter coefficient a, , relies on am m_ and we stipulate 
am,Q * ~1, when k>m, an , = %, the equation (11) is established for all m. 
(m+1) prediction error filter output power P,, is 


_ 


* tanm EUCE es ta} A Eevee y 





re — E C( Bas dar aBar)’ + (Bai— @eaBei)*) (12) 

in which 
Daw aber > Za. tom tbren (13a) 
Bai= Das saben o = Baw inn ebien (13b) 


. ¥ 


The (m+l) term prediction error filter's average output power minimization 
conditions dP,/da.,.-0, expressed by equation (12) have 


59 











> ¥ BaiB,. 
@asc —- (14) 
dy ‘8B: + Bo: 


From equations (11) and (13) we obtain the recursion formula of Bj, Bri 


B.. B.. lei Be ice iB le (15a) 
I S lese™ om 
BL. - B. le:cel ~Be lee 1B. lee (15b) 


The number groups B,, and Boi are made up of By_) and Bn-1,i multiplied by a 
simple linear operation, the initial value Bo; = Boy = ¢q and the value of 
Byj4 and Bj; is 


B:.=6 (j=1,2,°"",n) (16a) 


B).= Gi.: (i= 1,2,°-.n- 1 (16b) 
Substituting equation (11) in equation (7) we get the recurs: on formula of P, 


P.=P. i(l-az .) (17) 


Using equations (11), (14), (15), (16), and (17) we can directly solve the 
prediction error filter's coefficient ay, and the prediction error filter's 
minimum average output Pp directly from the raw data (), C25 «++s Sn» and 
estimate the ocean wave spectrum by equation (6). Represented as a circular 
frequency, equation (6) has 


P.At (18) 


the sea wave spectrum is generally expressed as a single-sided spectrum 
E(w) = 2S(w). Then the sea wave factors can be calculated through the 
spectral moment M, 


Mia | wo E(odo (19) 
the average wave height: H=, 2nM. (20) 
9 ° T ” Mo (21) 
the average period: eX J- d 
, M: 
the spectral width: ) j MoMe— i (22) 
MoM, 
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The length m(<n) of the prediction error filter influences the quality of the 
MEM estimate, determining the filter length m proper is an important issue in 
MEM estimates. We used the Akaike (1969a)°, 1969b9 and 197019 fixed Stage 
method and capacity ratio tixed stage method to determine the prediction 
error filter length for the different kinds of information of wind-wave, 
swell, and mixed wave. In this paper, the Akaike step determining method 
took the first very small value of the minimal final prediction error (EPE) 
and for centralized sampling we had 


EPE). em 
N—-(m: 1) P. (23) 
in which N is the sample capacity, m is the prediction error filter length. 
In the capacity ratio fixed stage method, we took 
m= N/20 (24) 


[I1l. Comparison of Spectral Estimates 


We used different information materials (wind-wave, mixed wave) observed by 
different wave measuring instruments (water pressure type, gravity type, 
ultrasonic type). We applied auto-correlation, period mapping, and MEM for 
spectrum estimating and in MEM aiso took the prediction error filter length mo 
determined by equations (23) and (24). 


1. Comparison of Methods for Spectrum Estimating 


Figures 1, 2, and 3 give the wind-wave materials observed by the water 
pressure wave measuring instruments, all used a sampling step length of 

\t = 0.375 second, and the spectrum curve of the various methods of esti- 
mating the spectrum where sample capacities were N = 2048, 1024, and 512, 
respectively. In the figures, COV represents the auto-correlation function 
method, FFT represents the period mapping method, MES represents MEM. In 

the tigures, in the conventional methods, k represents the free degrees, in 
MEM it represents prediction error filter length k, Akaike’s fixed stage 
method is determined by the self selected first extremely small value, 
therefore the value k is not clearly noted. In the figures, the solid line 
is the spectrum curve of Akaike'’s fixed stage method MEM spectrum estimate. 
From the three figures one can see that the MEM extreme value is higher than 
the extreme value given by the conventional spectrum estimation, and that the 
curve is smooth, indicating that the spectrum deviation is small. However, 
in Figure 1, when N = 2048, the extreme value of the Akaike fixed stage 
spectrum estimation is lower than the spectrum value estimated by the 
capacity ratio stage fixed at N/20 method. From the figures, it can be seen 
that MEM spectrum estimating is generally superior to the conventional method, 
the application of Akaike's method of determining m by the criterion of the 
first extremely small value is feasible in wind-wave information, because the 
two methods of determining m in MEM, when capacity N is large, the capacity 
determinant ratio method displays better spectrum extreme values, but the 
high frequency spectrum values are not as smooth as the Akaike fixed stage 
method. From Figures 1 to 3, when the imple capacity N shrinks to N/2 or 
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N/4, the spectrum values estimated by MEM have a spectrum curve similar to N 
indicated the superiority of MEM for short series spectrum estimating. 
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Figure 1. Spectrum Esti- Figure 2. Spectrum Figure 3. Spectrum Estimate 
mate (water pressure wave Estimate (water (water pressure wave measure- 
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N = 1024) Key: 
Key: 1. Legend 
l. Legend Key: 2. Akaike fixed stage 
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Figures 4, 5, and 6 give the wind-wave materials observed by ultrasonic 
wave measuring device. The step length for all samples was At = 0.75 
seconds, and the spectrum curve of the various methods of estimating the 
spectrum where sample capacities were N = 1024, 512, and 256, respectively. 
The solid line in the figure represents spectrum curve of the MEM estimate 
with capacity ratio stage fixed at N/20. Discussing the quality of the 
spectrum estimate from the angle of the spectrum's piandu [0252 1653 bias?] 
and deviation, we see that the extreme value of the spectrum estimated by 
the capacity fixed ratio method is the highest, but the fact that degree of 
freedom k = 20 in the high frequency spectrum value curve and the auto- 
correlation function method is rather uniform means that the capacity 

fixed ratio N/20 for estimating the spectrum values is rather good. The 
broken line in the figures represents the spectrum curve of the Akaike 
fixed stage MEM, the high frequency partial spectrum value deviation is 
small, and it is also superior to the spectrum curve of the estimate of 
degree of freedom k = 20 of the auto-correlation method. Therefore, Akaike 





tixed stage MEM estimated srectrum is better. When the sample capacity is 
reduced by doubling [chengbei 2052 0223], the MEM spectrum estimate give an 
excellent spectrum curve. 
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Figure 4. Spectrum Esti- Figure 5. Spectrum Figure 6. Spectrum Esti- 
mate (ultrasonic wave Estimate (ultrasonic mate (ultrasonic wave 
measurement record, Wave measuremen measurement record, 
N = 1024) record, N = 512) N = 256) 
Key: Key: Key: 

l. Legend l. Legend l. Legend 
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Figures 7, 8, and 9 give the mixed wave materials observed by the gravity 
wave measuring device. The spectrum curve of the various methods of esti- 
mating the spectrum where the step length for all samples was At = 0.75 
seconds, and where sample capacities were N = 1024, 512, and 256, respec- 
tively. The solid line in the figure representing the spectrum curve of the 
MEM estimate at a capacity ratio fixed stage of N/20, is superior to the 
spectrum curve of the conventional method and the Akaike fixed stage MEM. 
The MEM spectrum curve given by the Akaike first extremely small value dis- 
crimination method is a spectrum that tends to be smoothed, and the spectrum 
extreme value is lower than the conventional method, the corresponding 
frequency position difference for the spectrum extreme value is also 
greater. This indicates that the prediction error filter length m determined 








by the Akaike first extremely small discrimination method is suited to spectrum 
estimates of the pure Gaussian processes, and for mixed wave information is 
given a smoothed spectrun. 
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Figure 9. Spectrum Estimate (gravity wave measurement record, 
N = 256) 


Key: 
l. Legend 
2. Akaike fixed stage 


54 





2. Comparison of Spectrum Moment and Wave Factors 


Table 1 gives the spectrum moment Mg, M2, M4, spectrum width and average 

wave height, and average period of the wind-wave materials observed by the 
pressure wave measuring device. In the first column, MESA represents the 
calculated spectrum moment and wave factors by the Akaike fixed stage MEM 
estimated spectrum, MESB represents the capacity ratio N/20 fixed state MEM, 
COV is the auto-correlation function method, FFT is the period mapping method. 
From Table 1 it can be seen that sample capacity N = 2048, 1024, 512 correspond 
to the spectrum moments Mg, M2 and average wave height H, average period T, 
and spectrum width calculated by the various methods are all fairly close, 
only the relative changes in My, are larger. However, for actual application, 
the above described parameters are sufficiently precise, therefore the 
spectrum moments and wave factors estimated by the various methods have 
practical significance. The wind-wave materials measured by the ultrasonic 
wave measuring device also produced similar results, and they will not be 
tabulated and explained here. 


Table 1. Comparison of Spectral Moment and Wave Factors (wind-wave materials) 






































(ly # (25RN Mo M: Ma, e H T 
MESA 2048 2° 298 2.017 16.172 0 944 3.800 6.71 
MESB 2048 2.377 2.081 16.280 0.942 3.861 6.71 
COV 2048 2.444 2.101 13.905 0.933 3.918 6.78 
FFT 2048 2.407 2,094 15,404 0.939 3.889 6.74 
MESA 1024 1.890 1.750 10,267 0.937 3 450 6.53 
MESB 1024 2.070 1,878 13.308 0.934 3,607 6.60 
COV 1024 2.040 1,810 11,043 0.924 3.581 6.67 
FFT 1024 1.998 1,808 13,044 0.935 3.540 6.61 
MESA 512 2 060 R40 11.410 0.937 3.60 6.65 
MESB 512 2.004 1.808 13.306 0.987 3.55 6.62 
COV 512 2.020 1.750 10, 460 0.924 3.56 6.75 
FFT 512 1.919 1,712 12,589 0.937 3.47 6.65 

Key: 
1. Method 


2. Capacity N 


Table 2 gives the mixed wave materials measured by the gravity wave measuring 
device, and applies the spectrum moment, spectrum width and wave factors 
calculated for the estimated spectrum by the various methods. From Table 2 
it can be seen that the statistical values found for the spectra estimated 

by the various methods are all very ciose, and have practical significance. 
The statistical values (spectrum moment and wave factors) of Tables 1 and 2 
explain that whether wind-wave or mixed wave, although there are great 
differences in the spectrum curves estimated by the various methods, the 
statistical values are very close and can be used for general engineering. 


65 





Table 2. Comparison of Spectral Moment and Wind Factors (mixed waves) 

















(1) . (2)#'N Me M; M, e€ H T 
MESA 024 0.321 9.45 912 0.818 1.419 5 30 
MESB 1024 320 > 449 1 901 0.813 1.418 5 30 
CoV 102 ) 312 0.445 1 885 0 817 1.412 5.30 
FET 1024 ) 31 0.434 790 0.813 1.396 5 31 
MESA 326 45: 1 898 0 820 1.453 5.39 
MESB 512 335 0.456 1 887 0.819 1 451 5.39 
COV 51 329 448 845 0 81% 1 437 5 3% 
FFT 512 9.301 0.421 742 0.914 1.376 § 31 
MESA 256 0 385 . 566 ° 855 0 843 1.555 § 19 
MESB 56 22 564 ° 865 0 843 1.555 5.19 
COV 256 0.370 ) 547 ° 752 0.841 §25 § 17 
FFT 256 0.344 0.515 2.532 0.841 1.471 5,19 

Key: 


1. Method 
2. Capacity N 


IV. Conclusion 


1. For wind-wave materials measured by water pressure wave measuring 
devices, the spectrum curve estimated using the Akaike first very small 
value fixed stage MEM is superior to the spectrum curve estimated by the 
conventional method and the capacity ratio N/20 fixed stage MEM. 


2. For the wind-wave materials measured by the ultrasonic wave measuring 


device, when the sample capacity N > 512, the spectrum values estimated 
using the capacity ratio N/20 fixed stage MEM are better, and it is superior 
to the conventional method and the Akaike fixed stage MEM, when N = 256, the 
spectrum values estimated by the two fixed stage MEM are close. This indi- 
cates that when the ultrasonic wave measuring records which directly reflect 
the wind-wave information are used to estimate the spectrum, the smoother 
spectrum produced by the Akaike fixed stage MEM has a resolution lower than 
the capacity ratio N/20 fixed stage MEM and is superior to conventional 
method of estimating spectra. For convenience of use, we can use the 
Akaike fixed stage MEM. 


3, For mixed waves, the spectrum estimated by the Akaike fixed stage MEM 
vives a smoothed spectrum, resolution is lower than or equal to the conven- 
tional spectrum estimate, but the capacity ratio N/20 fixed state MEM esti- 
mated spectrum displays better resolution and is superior to spectrum 
values estimated by conventional methods. Therefore, when the Akaike fixed 
stage MEM is used to estimate the spectra of mixed wave materials, it is 
still necessary to improve the discrimination standards for determining 

the prediction error filter length. Otherwise, carrying out test 


66 





calculations using the capacity ratio fixed stage method is like trials 

using the conventional method with different degrees of freedom to discrimi- 
nate the optimum quality spectrum estimate. 

4. Although in using the various methods to carry out spectrum estimates 

on wind-wave (6 groups) and mixed wave (4 groups) wave surface records 
measured by different devices there are differences in the spectrum 
structures and difference deviations and biases are revealed, the statistical 
parameters--spectrum moment, spectrum width, average wave height, average 
period--through calculations of the estimated spectrum are all fairly close 
and all have practical significance. 


5. For short sample capacities, the Akaike fixed stage MEM estimated 
spectrum exhibited that it is even more superior to the conventional 
spectrum estimating method. 
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PHYSICAL SCIENCES 


PHASE SPEED OF WIND-WAVE COMPONENT STUDIED 


Qingdao SHANDONG HAIYANG XUEYUAN XUEBAO [JOURNAL OF SHANDONG COLLEGE OF 
OCEANOLOGY] in Chinese Vol 15 No 3, 15 Sep 85 pp 33-39 


[Article by Xu Delun [1776 1795 0243], Department of Physical Oceanology and 
Marine Meteorology: "Measurements of Phase Speed of Wind-Wave Component and 
Interpretation of Their Result"; article received 30 June 1984] 


[Text] English Abstract: Measurements of phase speed of wind-wave component 
were made by cross-spectrum method, separation distance between the two probes 
for measuring surface elevations being 2 cm. 


The measured data indicate that the relationship between phase speed and 
frequency varies with wind speed and water depth. For one extreme condi- 
tion, low wind and deep water, it obeys basically the classical linear 
dispersion relation. For the other extreme condition, high wind and 
shallow water, the components are non-dispersive, all traveling with the 
same speed as the component with spectrum-peak frequency. For medium condi- 
tions, the dispersion relations have some forms between the two extreme 
cases. 


The measured result can be interpreted on the basis of a refined spectra! 
model, which consists of the contributions of both dispersive and non- 
dispersive wave with the effect of drift currem included. At low wind and 
in deep water, the wave field is of low nonlinearity: free waves are domi- 
nant and, hence, the wave components are dispersive and obey the classical 
dispersion relation. At high wind and in shallow water, the wave field has 
high nonlinearity. In this case, forced waves dominate over free waves and 
the wave components appear to be non-dispersive. 


lL. Introduction 


Linear wave spectrum theory holds that wind-waves are the result of the 
linear superposition of simple waves of many frequencies, unequal oscilla- 
tion and confused phase. According to this theory, the phase velocity of 
wind-waves should obey linear wave dispersion relation c = g/n (n is circu- 
lar frequency, g is gravitational acceleration). Does this theory really 
conform to a large degree to the actual wind-wave conditions? The answer 
relies on accurate measurements of wind-wave phase velocity. Ramamonjiarisoa 
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(1974)1 of the Marseille Laboratory in France used the cross-spectrum method 
to measure and obtain data. This data indicated that only the wave component 
of the wind-wave whose frequency was near the spectrum peak frequency obeyed 
the linear dispersion relation, wave components with higher frequencies were 
non-dispersive, and both were propagated at almost identical phase velocities. 


The measurement data of Ramamonjiarisoa (1974) attracted wide attention. 

Some researchers (Huang, 1981; Dudis, 1981, et al.) expressed doubts about 
the measurement method of this data. Through a strict theoretical analysis, 
Huang (1981) pointed out that the distance between the two probes which 
Ramamonjiarisoa (1974) used in measuring wave surface height was too great 
(10 cm) and thus he could not measure the phase velocity of wave components 
at higher frequencies. In addition, Ramamonjiarisoa's measurement data was 
only obtained at medium velocity, and lacks data at low wind speed and higher 
wind speeds. 


This paper first of all reports the laboratory measurement results. These 
data were still obtained using the cross-spectrum method, but the distance 
between the probes was reduced to 2 cm, and the results include various 
situations such as from low wind speed (3 m/sec) to higher wind speeds 

(14 m/sec) and shallow water. Finally, there is a theoretical interpretation 
of the measured results. 


Il. Cross-Spectrum Method 


The definition of cross-spectrum ¥(r,n) is 


Wren) PF D (ken) aa dR (1) 


* * 


in which, $(R,n) is the wave number-frequency spectrum, k is the vector wave 
number, n is the circular frequency, r is the vector distance between the two 
measuring probes. 


The phase spectrum function can be found from the cross-spectrum 


QC r.n) : (2) 


C,(r.n) 


Xr om) - tan 


in which Cg and Q are the real parts (equal phase spectrum) and false parts 
(out-of-phase spectrum), of the cross-spectrum ¥(r,n), 1.e., 


a —_ 


Col r.n)= |u.cos(k. r)oCk. nddk 
. 
and 


OCr, n)=| sini k.t )O¢ k,n Ddk 


I0 








it the wind-wave is propagated only in the single direction r, then the phase 
velocity spectrum of the wind-wave component can be found from the wei [0143 
positional?] phase spectrum 


(n)= 22 
6(r,n) (3) 


- , ; _ , _ | - , 
ihis equation is the one used by Yefimov et al. (1972)- and Ramamonjiarisoa 
et al. (1979)3 to calculate the wind-wave component phase velocity. 


Atter stricter theoretical deductions, Huang (1981}* stated: The suitable 
conditions tor equation (3) only limited to wave field is unidirectional, 
and the distance between the two probes must be small enough to satisfy the 
condition 


Rerl «i (4) 
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lll. Experiment and Results 


Experiments were carried out in a wind-wave-flow experimental trough at the 
t\iversity of Delaware in the United States. The trough length was 44 meters, 
| meter wide and the water depth at the time of the experiment was 0.75 meters. 

[wo capacitance wave measuring devices were used for measurement of the wave 
surtace height. The distance between the two probes was 2 centimeters. Four 
groups of measurement data were obtained in all. Three of the groups were 
taken in a wind area of 30 meters, at wind speeds of 3, 8, and 14 m/sec, 
respectively; the other group was taken in shallow water on a slope at the 

end of the trough, wind speed was 12 m/sec. The recording time of each group 
f data was 304 seconds. The sectional direct method was used for estimating 
the frequency spectrum, each 304 second record was broken up into 38 sections. 
fhe sampling rate was 64 Hz, there were 19456 data points in each group, and 
912 data points in each section. The frequency (f = —) resolution of the 


cross-spectrum used for estimate results was 0.125 Hz, and its degree of 


I reedom was 76. 


Since the water trough was narrow, we could more or less accept that the wave 
field measured was unidirectional. Probe distances of 2 cm make the wave 


aponents below a frequency of 5 Hz satisfy condition (4) very well. This 
means that the results of the measurement are effective for wave components 
f below a frequency of 5 Hz. The phase velocity relative frequency f points 


ilculated according to equations (2) and (3) from the wave surface records 
ire plotted in Figures 1, 2, and 3. For ease of comparison, two separate 
irves are also plotted in these figures. The solid Line is the theoretica! 
irve ileculated from the linear wave dispersion relation (including the 


’ 


influence of surface tension T). 
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Figure 1. Relationship of Frequency and Phase Velocity of Wind-Wave Component 
at Low Wind Speed (Wind speed is 3 m/sec, water depth is 75 cm) 


Key: 1. Phase velocity (cm/sec) 
2. Frequency (Hz) 
3. Standardized spectrum density 
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Figure 2. Relationship of Frequency and Phase Velocity of Wind-Wave Component 
at Medium Wind Speed (Wind speed is 8 m/sec, water depth is 75 cm) 


Key: 
1. Phase velocity (cm/sec) 
2. Frequency (Hz) 
3. Standardized spectrum density 
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Figure 3. Relationship of Frequency and Phase Velocity of Wind-Wave Component 
at Higher Wind Speed (Wind speed is 14 m/sec, water depth is 75 cm) 
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o =| BA O- 0’, _2xT (5) 





(in which op and 9" are water and air density, respectively and A is wave 
length); the broken line was obtained from the right side of equation (5) 
plus the phase velocity c, which increases because of wind drift, i.e., 
from the equation 





,_ Ip ga e-0’ 2xT 
© "Vox o+n Apre) ° (6) 





In the laboratory, the drift produced by the wind had a very clear impact on 
wind-wave velocity. We do not yet have a theoretical method to calculate 
this impact, nor do we have actual measurement results which can be used. 
For ease of comparison, we use approximations [7], we view the differences 
at the spectral peak frequencies of velocity measured by the cross-spectrum 
method and the velocity provided by the linear dispersion relation as wave 
component phase velocity of value cy which increases because of wind drift. 


From Figures 1-3 it can be clearly seen that at different wind speeds, the 
tendency of wind-wave component phase velocity which was measured changes 
with frequency are slightly different. At low wind speeds (3 m/sec), the 
measurement data of frequencies below 5 Hz fits fairly well with the 
theoretical curve. At medium wind speeds (8 m/sec),. beginning with a 
frequency of 2.3 Hz, the measurement data systematically deviates from the 
theoretical curve. This is generally the same as the results measured by 
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o( k,n) = @x(n) &k-—W—) + oo(m) {u-— Fr) (7) 


in which >y and Pp represent, respectively, the contributions of non- 
spersive and dispersive waves to the frequency spectrun, Cy is the phase 
velocity of the non-dispersive waves, c,, is the phase velocity of the wave 
component increased due to wind drift, and 6 represents the Dirac delta 
tunction. 


sudstituting equation (8) in equation (2) we get 





rn . 
Massie. +s 
Cu +c +c 
tan8(n) = a —_— (8) 
u(n)cos—"_ + cos 2 — 
Cw rec, g*c.n 
in whicn 
(mn) = Oxi) (9) 
on(n) 
is the ratio of the contributions to the frequency spectrum of non- 


sispersive waves and dispersive waves. Ina research report, Xu et al. 
(1984) solved u(n) as a function of effective wave surface slope ¢ = 


(in which, (€*)'”* is the root-mean-square surface height, \g is the 
wave length of a wave component which has the spectrum peak frequency) and 


the effective depth (n is the water depth). For deep water and low wind 
speed situation, u(n)<<l, thus, equation (8) is simplified to 





rn | 
vin) - 10) 
g*c.n 
Substituting equation (10) in equation (3) we get 
cin) Ss * ¢, (11) 


n 


Clearly, under low speed conditions, except for the added term c,, produced 
by wind drift, the phase velocity of the wave component fits to the linear 
dispersion relation equation. For shallow water and high wind speed, u(n)>>l, 


thus, equation (8) is simplified to 


Fn (13) 
cx +c [sic] 


M(n) 
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-HINESE LANGUAGE MICROCOMPUTER DESCRIBED 


Seijing DILANZI KEXUE JISHU [ZLECTRONICS SCIENCE AND TECHNOLOGY] in Chinese 
1 16 No 5, 10 May 86 pp 10-11 


(Article by Chen Shiming [7115 0013 6900], Academy of Sciences, Guangdong: 
The Characteristics and Use of the LX-80A Chinese Language Computer” j 


Text) The LX-80A associative Chinese character graphics microcomputer system 

(hereafter, Chinese language computer) was successfully developed by the 

Ins? te of Computing Technology of the Chinese Academy of Sciences, and ‘s 
irrently in production by units such as the experimental factory of the 
uuangdong Academy of Sciences. Its characteristics and methods of use are as 


follows 
Multiple Chinese Character Input Methods 
Chinese character input for the LX-80A Chinese language computer is easy to 


learn and is broadly adaptable. 


he “interactive” Chinese character input method. 
rn method is ded according to the stroke shapes of Chinese characters, so 
it is also called “stroke coded.” There are 40,000-50,000 Chinese characters, 
more than 6,000 of which are in common use, and stroke coding can classify the 
stroke components of the character radical into 26 basic shapes. There are 
6 basic strokes (~|, » £7) , 7 common radicals (4, ®. 7. #. A. #. ), 


12 stroke combinations, and 1 independent character-stroke shape. ihe code of 
each stroke uses the keyboard position of an English letter, and the name of 

each stroke uses the Chinese character command that may be seen at the left of 
Figure 1. For example, the stroke shape ‘tt’ is called “grass and the corres- 


ponding key is the letter Q. 








| 


=_—r 1 


Wet The /iig@ “lla All tlhe Clian &I 
cM aj 7 gi of @ oP a Aw(* vox on 4 








/ ————.}« 


it se! sit TT if . 
| a 
Sz tr zcq@v tetn te 





1 on a 2 , aad Ba<s ] Vv ‘\ 4 : mm < cak 1 Se ee, a eee | nal a oa 
ine input method used Dy the LA-SUA System 1S Calied Visual Keyboard input. 


‘ +} ha ~ £ 1 _ — . : _ < | So (4 lin nramnt 
At tne pottom or tne screen tf the display device there is a three line prompt 
irea, where after power up, the 26 ctroke shapes in Figure 1 and their 


7 smravasla , Bnolicl 1 7 Is 1 , ee a 1- , r a . 
-Orresponding Engiisn ietters are ail displaved in the rompt lines. io enter! 


a particular Stroke, just press tne corresponding tetter key. 10 enter Cne 
" : _ 1 7 ; 7 - _ tir oo 
Prass radicai, press VW. 
7 
US y Stroke snape OdIL!I to enter Uhines naracteyr is iargely the same as 


rroceeds trom left to right rrom top t DLULCLUM trom the outside to the 
; i ’ 

inate vee The stroke shape coding method of the LX-80A svsten . 

inside, etc. he stroke Shape coding metnod oOIf the LA-dSUA system use 


& 2 =. ‘ ~ ] 7 | +} . j . F . - | . + } 
tne first Stroke ind the second stroke, Those characters that use on] the 
irst stroke may then b Tound, Dut the ma rity of characters will need the 


S 
econd stroke shape bDetftore they can be tround,. ror some words that are not 
—- ¢ . nr, | . ” ; . } ; : 
ent words, as for example a » after the first str 
ty 9 ? 


te = , 24 3 
IViISiIDIe, cal led independ 
ipe press the L Key (C 


; 4 . ow 4 
he independent word key), which will then locate 
the word, [The distinction between the first stroke shape and the second is 
11.adicated DY the letter appearing if the ieftt in the prompt area, we Wilil NOW 


1} h ' ear vf | c ha . c¢ Mots. ‘ , Cartas " Te 
use the entry ot the )») Character oO] Chinese Academ I lences fe] BF en 


5 | 
. . ++ + +} , +} + ¢» ' 
A 4 — my at CT. LUSLCITACE eiba So See | . at Il! >t ne nme a7 | 
> > : + - ? ’ 
) il Cers is S 1] the writ l .a ] ivure . 
? —7 «+,)} , = de 
rer a. an —_—_—* — 
se. 7— Le = \ 
} _ — 
‘ —_— ‘4 ‘ 
} 1Y 
A Ai - 
. ' } ' ' te I ‘ 
: > ‘4 mer* + = + . . > + rm > + ; ; 
le Llirst stroke oO! zn i. is] if » Key. tC this time, 


haracters havin : ’ eP.r tirst stroke will appear in the prompt area 


screen (the 26 characiu.s In the prompt i1rea Nave dittering stroKe unt 


] mar ¢ y ‘ + Lx } +} — = M y } Tr e } ~ ha - ar » ’ >} \ ’ 
iot necessarily be the same eacn time), and when the character ZNONE 
. ’ ‘ tae ft . tr, 
ap i t he sition of the lL : ress lower cas M" to enter the ' 
. ' ! F . —_ ' ' ' ’ : 
narecte 2non 1 the text. l , Di e ¢ iracte! yuU¢ ress N,. 
' ' ' ' : ' ' ' ' 
then enter tl lower Cast ; r t haracter ke’ press 'F, for ut 
' ' ' ’ ' F . . - ’ ’ ‘ 
pre y ind fOr uan re ; ter hich 1 must also pres ) for 
t rtrn ? y ' ' - ’ > , 
er G L ] Ke Snape ( re yu i ear L rie SCTe? ° 
) t ] f } ) ; ‘ ] " + } y ‘ I hi 
»¢ J f i ISI1NZ Z0 Stroke snape [ l ] j ilne ( ° i) 
T ¢ i I ¢ ite [ ré ; Te tt nt if Lon t i 1cit i j { I nt n ait lL] \ nil it ¢ 
! Le! ° I le, t} 7 } CSE { pDOtN *', maa ) 
” - = : P ; ? 
} < CULCe, r tn { i} I these nay ters all have ri nape j 
. - , | ¢ . . y i; y / 
ising Stroke e codil enter CUne cter ! l r the qgica 








2. The romanization entry method. 


Similarly to the “interactive” entry method, use the romanization key for the 
LX-80A system to accomplish the task of character entry by pressing the keys 
that represent the initials and finals and the corresponding material in the 
prompt area of the screen. 


3. The “arbitrary selection coding" entry method. 


As of this moment, there are hundreds of different Chinese character input 
methods, and the design of the LX-80A has taken into consideration coding 
schemes with which the user may already be familiar, as for example the tele- 
graphic code, the Cang Jie code, or the four-corner code, and the user may 
select one among them as an input method for the LX-80A system. After entering 
the telegraphic code, press the “arbitrary selection" function key and the 
computer will be placed in the “arbitrary selection coding" state (i.e., 
telegraphic code). At this time, characters may be entered with telegraphic 
codes, 


The LX-80A system provides users with software for “automatic generation of an 
arbitrary selection coded dictionary’ that will allow the user to use the 

coding method with which he is most familiar or to change to a new coding scheme 
designed by the user. 


4. The “international index code" entry method. 


Chinese character international index codes are four-digit hexadecimal codes, 
used, for example, by sending a "30AE" to represent the character for "to 
love."' The transmission steps are as follow: press the “international index 


code" function key and the computer will then be placed in the international 
code status. Then, enter a "30AE" and press the space bar. The computer will 
display the character for "to love." If more than 4 digits are entered, then 
the code will automatically be discarded, only the last four codes being 
accepted. With international coding, Japanese and Russian may also be entered. 


The four entry methods just described are controlled by corresponding function 
keys, and methods may be changed at any time for greater convenience. 


Independent Chinese Character Associative Functions 


There is an associative dictionary in the LX-80A system that contains tens of 
thousands of common words and phrases. As each Chinese character is entered, 
the computer guesses the next possible word based on the dictionary, and then 
displays this word on the screen prompt area. The operator may then directly 
hoose this word. A character will bring up a group of characters (26), which 
is very like the human associative function. In this way, and by interlocking 
associations, a large number of words, phrases, and even whole sentences can 
be generated. The partia] associative structure for the character ‘zhong’ in 


? 


Figure 3 is shown with arrows indicating the directions. 
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Figure 3 


Use of the associative function is very easy and may be divided into the 
associations with "characters" and the associations with "words." 


1. Associations with "characters." 


This associative function is controlled by function key F7. Pressing this key 
will present the word “association” in the prompt area, indicating that the 
computer has entered the associative state. If the character "zhong" is 
entered, 26 characters like "hua" and "guo" will appear in the prompt area. 
After choosing "guo," 26 words like "qi," "jia," and "ke" will appear in the 
prompt area. In this way, entry of the five characters "Zhongguo Kexue Yuan" 
("Chinese Academy of Sciences''] can be by the associative method and need not 
be entered by the character at a time method described earlier. 


2. Associations with "words." 


This is controlled by function key F9, and if after entering the character 
"zhong' F9 is pressed again, many words and phrases associated with the 
character 'zhong' will appear in the prompt area, as for example, "leading 
comrades of the Central Committee," ''Chinese Academy of Sciences," etc., and 
when there are many associated words and phrases, they may be displayed by 
repeated pressing of F9. If a particular phrase is to be entered, the key 
corresponding to the final character in the associated phrase may be pressed 
and the entire word or phrase will be entered into the computer. The longer 
the sentences are, the higher the efficien 


There is associative function support for each of the character input methods 
7 


described above, which allows character input speed to reach an average of 
1.5 key presses per character. 





Command operation of the LX-80A system itself is now by Chinese characters. 
wenty-six Chinese character commands have been burned onto a 64 kb EPROM, 

ind upon system reset, the 26 character commands are all displayed in the prompt 
irea by their functions. As shown in Figure 4, if the user wishes to see the 


‘ry 1 


tile directory of the current disk, he may press Q, and the screen will then 
immediately display the file directory and will also tell the user the size of 


each file. If one wishes to know which erations commands are represented b 


. 
Cu 
f 

r 


the English letter keys, press ‘W screen will display the meaning 


each command. 


The LX-80A is a multi-function, stand-alone Chinese language computer systen 
having 576 kb of RAM and two flopp: isk drives, with a 24 pin printer. It 
can be used for document processing, storage and display, and for printing. 

it is equipped with a standard interface (RS-232) and may be connected with a 
variety of host computers, which makes the LX-80A a powerful Chinese languags 


terminal. 
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INVESTIGATION OF UNDULATOR FOR HEFEI 800 MEV ELECTRON STORAGE RING 


Beijing GAONENG WULI YU HE WULI [PHYSICA ENERGIAE FORTIS ET PHYSICA NUCLEARIS] 


Article by Chen Zihcheng [7115 1807 6134] and Liu Jian-yang [0491 1696 2254], 
China University of Science and Technology, and Xu Jianming [1776 1696 6900], 


ryt ° . o . tt, : 
Physics, ULninese Academy of Sciences: The Investi- 


Institute of High Energy 
for the Hefei 800 MeV Electron Storage Ring"; paper re- 
, 
j 


r 
gation of an Undulator 
ceived on 11 March 198 
[Text] Abstract: A permanent magnet undulator for the Hefei 800 MeV electron 
storage ring is now under construction. This paper introduces undulator radi- 
ation characteristics and magnet design considerations and reports on the de- 

Sign, development, and magnetic field measurement results of a 10 period hybrid 


“4 ¢ ] 
unduiator model. 


An undulator is a periodic magnetic device which, when inserted into the elec- 
tron ring's linear section, can supply quasi-monochromatic synchronous radia- 
light of greater brightness than the radiation produced by a wanzhuan 

37 6567 turning?] magnet, and thus becomes an important insertion component 
storage rings. Its technical reliability has been verified by successful 
experiments by many laboratories abroad.* The undulator has the following three 
primary features: 1) The special structure of the undulator magnet makes the 
ring's electron beam not produce net deflection or shifting and does not have 
1 serious influence on the electron beam movement; 2) In the wanzhuan magnet, 


| 


the electron beam radiates a ntinuous and smooth light spectrum, but in the 


> 
) 


() 


undulator, the interference effect of the radiation suppresses the radiation 
light in some wavelengths, thus in s1 ific wavelengths it is intensified, and 
the radiation is concentrated within a narrow wavelength interval, and the 
light spectrum obtained is several narrow quasi-monochronomatic light peaks; 

3) TI radiation focuses within a very small space solid angle along the axis, 
and has a high spatial brightn 3 and llimation, which reduces the heat given 
off by the beam ray mponents, and radiation power can be more effectively 


When high magnetic energy rare earth metal permanent magnet material is used 
for the undulator magnet, an excitation current is not required, there is no 
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Figure 3. Horizontal Directivity Chart Figure 4, Vertical Directivity Chart 
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Figure 5. Relation Curves for Standing Figure 6. Diagram for Expanding 
Wave Ratios and Frequencies Bandwidth 


2. Testing frequency characteristics. 


For these tests we selected channels 6 through 12 in the VHF barid. The 
frequency characteristics of the half-wave square--loop antenna were indicated 
through tests of the standing wave ratio. The antenna input impedance was 
300 ohms, and the antenna input terminal was connected through a broad band 
impedance transformer and 75 ohm coaxial cable. 
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Using a square-loop antenna made of copper wire of differing diameters, the 
relational curve of the standing wave ration S and the frequency ~ is as shown 
in Figure 5. 


Curve 1 shows the design frequency F ,=195 MHz of the square-loop antenna, where 


side A =A /2=770 m and the diameter of the antenna is p = 4.4 mn. 


Curve 2 shows the design frequency f 97195 MHz of the square-loop antenna, 


where side 4 =,A/2=770 mm and the diameter of the antenna is p =1.5 mm. 


Curve 3 shows the design frequency Ff ,=203 MHz of the square-loop antenna, 
where side & =) /2=740 mm and the diameter of the antenna is $ =8 mm. 


Ill. Discussion 


l. There is a definite difference between the tested horizontal directivity 
charts and the horizontal directivity charts by theoretical analysis. -The 
direction with the greatest dip just happens to be the direction where the 
signal is aligned with the side that meets the feeder line. This shows that 
for the side that connects with the feeder line, reception efficiency is not 
as good as for the sides without the feeder line. In two tiered applications 
of square-loop antennas, this condition should be in mind. 


2. It can be seen from Figure 5 that the higher the frequency of the square-loop 
antenna application, the thicker should be the diameter of the material used. 

The curve 2 in Figure 5 shows that when the diameter P| #1.5 mm, the low point 
(179-187 MHz) standing wave ratio of the test frequency is smallest. Curve l 
then shows that when the diameter 3 =4 mm, square-loop antennas of equal 
dimensions will have the least standing wave ratio at the high end (195-211 MHz). 
Curve 3 shows that when the antenna diameter used is g =8 mm, the practical 
frequency is even higher. 


3. To further develop the working frequencies of the square-loop antenna, the 
method shown in Figure 6 may be used for adding shorter reception antennas. 
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Electronics 


THEORETICAL, EXPERIMENTAL STUDY OF CAVITIES IN TWT 


Beijing DIANZI XUEBAO [ACTA ELECTRONICA SINICA] in Chinese Vol 14 No l, 
Jan 86 pp 7-15 


[English abstract of article by Zhang Zhaohong [1728 2507 3163] and 
Wu Hongshi [0702 7703 6684] of the Beijing Vacuum Electron Devices Research 
Institute] 


[Text] An analysis of coupled cavities in TWIT is presented. The analysis 
employs field equivalence principle in conjunction with Green's function 
techniques. The resulting operator equations are solved by the moment method. 
The characteristics of the coupled cavities are expressed in terms of general- 
ized admittance matrices which are functions both of frequency and dimensions 
of the system. The results of calculation performed with the aid of an elec- 
tronic computer agree quite well with those of cold test measurements. (Paper 
received November 1984, finalized July 1985.) 
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DYLMOS INTEGRATED CIRCUITS: HYBRID INTEGRATED CIRCUIT WITHOUT ISOLATION 


Beijing DIANZI XUEBAO [ACTA ELECTRONICA SINICA] in Chinese Vol 14 No l, 


Jan 86 pp 16-20 


{English abstract of article by Wang Shoujue [3769 1349 6030], Xia Yongwei 
{1115 3057 0251], Kong Lingkun [1313 0109 0961], Li Yuanjing [2621 6678 6975], 
and Ho Naimin [0149 0035 2494] of the Semiconductor Institute, Chinese Academy 


of Sciences, Beijing] 


[Text] A new hybrid integrated circuit 
developed. By increasing and adjusting 
tors and bipolar transistors with large 
AND-OR gate in DYL are developed on the 
characteristics among the linear AND-OR 


without isolation, DYLMOSIC, has been 
the technology, P channel MOS transis- 
8 and small 8 value for the linear 

same wafer, based on the self-isolation 
gates, the main logic unit in DYL, and 


MOSIC. Using the hybrid integrated technology, an inverter unit is constructed 
by bipolar and MOS transistors on the same N type silicon wafer. The benefits 
of DYLMOS hybrid integrated circuits are: simple and compatible technology with 
the linear AND-OR gate in DYL, high input and low output impedance and compati- 
ble level with DYL and TTL circuits. (Paper received March 1985, finalized 


July 1985.) 
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INFLUENCE OF MISALIGNMENT OF ELECTRON LENSES ON IMAGE-QUALITY 


Beijing DIANZI XUEBAO [ACTA ELECTRONICA SINICA] in Chinese Vol 14 No l, 
Jan 86 pp 21-29 


[English abstract of article by Tong Linsu [4547 2651 1121] and Wu Changjin 
[0702 1603 3160] of the Nanjing Institute of Technology, Jiangsu] 


[Text] The misalignment of the rotationally symmetrical electron optical 
system caused by manufacturing and assembling tolerance is an important limi- 
tation factor for further improving the image quality. The finite-difference 
method is first used to calculate the three-dimension field to analyze the 
problem of the misaligned electron lense. The influence of the misalignment 
of equidiameter cylinder bipotential immersion electron lense on image quality 
is given. Experimental tubes are also manufactured. It is shown that the 
theoretical computations are in good agreement with experimental values. 
(Paper received September 1984, finalized March 1985.) 
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NEW METHOD OF MEASURING TWO-PORT DEVICES BY SIX-PORT REFLECTOMETER 


Beijing DIANZI XUEBAO [ACTA ELECTRONICA SINICA] in Chinese Vol 14 No 1, 
Jan 86 pp 44-50 


[English abstract of article by Sun Jia [1327 4471] of the Chengdu Institute 
of Radio Engineering] 


[Text] A new circuit is given for measuring scattering parameters of any two- 
port microwave network by means of a single six-port reflectometer. Based on 
the wave superposition theory and the theory of equivalent circuit, a new 
method of calibration and measurement is described. (Paper received July 1984, 
finalized July 1985.) 
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CONSTRAINED OPTIMIZATION DESIGN OF ELECTRON OPTICAL SYSTEM USING COMPLEX 
METHOD 


Beijing DIANZI XUEBAO [ACTA ELECTRONICA SINICA] in Chinese Vol 14 No l, 
Jan 86 pp 51-57 


[English abstract of article by Gu Changxin [7357 2490 9515] and Shan Liying 
(0830 5461 5391] of the Modern Physics Institute, Fudan University, Shanghai] 


[Text] The constrained optimization method-complex method described in this 
paper is better than the simplex method for design of electron optical systems. 
This method is beneficial to making optimal automatic search and can ensure 
that the computed results will satisfy the design requirement. In the complex 
method as well as simplex method, it is not necessary to know the explicit 
functional relation between the objective function and the searching para- 
meters. Comparing with other optimization methods, the complex method and 
simplex method have significant advantage in the optimization design of elec- 
tron optical systems. (Paper received May 1984, finalized February 1985.) 
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ON-LINE FAULT DIAGNOSIS OF ANALOG IC 


Beijing DIANZI XUEBAO [ACTA ELECTRONICA SINICA] in Chinese Vol 14 No l, 
Jan 86 pp 65-72 


[English abstract of article by Zhao Guonan [6392 0948 0589], Xie Shuang [6200 
7208], and Guo Yushun [6753 5940 7311] of the Hangzhou Institute of Electronics 
Engineering, Zhejiang] 


[Text] The following new ideas are presented for use in on-line fault diagno- 
sis of analog integrated circuits. 


1. The number of fault elements in an analog circuit can be evaluated from 
the measured data upon the limited number of accessible terminals. 


2. The dependent relation between any two multiple-fault sets can be classi- 
fied in two kinds: the mutual dependence and the single dependence. The 
routine of on-line fault diagnosis can be accelerated with the use of mutual 
dependence. 


3. The idea of effective measuring ports is helpful to the topological 
judgment of diagnosability of certain multiple-fault set besides the examina- 
tion of the existence of loops and cutsets among the fault elements. 


4. The uniqueness of the values of certain existing fault elements calcu- 
lated from the measured data upon the different measuring ports can be used 
to distinguish the actual fault from the guessed fault. 


Program based on the algorithm derived from the above ideas has been proved 
with success. (Paper received October 1984, finalized April 1985.) 
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STATISTICAL COMPUTATION OF COMPRESSION RATIO FOR ZERO-ORDER PREDICTOR WITH 
FLOATING-APERTURE ALGORITHM 


Beijing DIANZI XUEBAO [ACTA ELECTRONICA SINICA] in Chinese Vol 14 No l, 
Jan 86 pp 72-79 


[English abstract of article by Guo Guirong [6753 2710 5554] and Shen Qiang 
[3947 1730] of the Changsha Institute of Technology, Hunan] 


[Text] The problem of statistical computation of compression ratio for zero- 
order predictor with floating-aperture algorithm is considered. Markov Chain 
is proposed as the model of input process. By this, analytical expressions 
are derived for the mean and mean-square value of compression ratio, andthen 
a simple and convenient algorithm is developed for the computation of these 
factors. Finally, several numerical results are given to demonstrate the 
effectiveness of this approach. (Paper received July 1983, finalized July 
1985.) 
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THEORETICAL ANALYSIS ON INTERMDDULATION NOISE IN ANALOG MICRC'JAVE TRUNKS 


Beijing DIANZI XUEBAO [ACTA ELECTRONICA SINICA] in Chinese Vol 14 No l, 
Jan 86 pp 80-86 


{English abstract of article by Yang Yunnian [2799 6663 1628] of the Fourth 
Research Institute of the Ministry of Posts and Telecommunications, Xian] 


[Text] A theoretical formula to explain the intermodulation noise, which is 
produced by general transmission characteristics with higher terms in an 
analog microwave trunk, is obtained by use of the Bessel function expression. 
The usual equations for calculating the intermodulation noise in engineering 
practice are special cases of the presented formula. To analyze a general 
transaission characteristic with higher terms in microwave trunk, the least- 
squared principle is recommended and satisfactory results are obtained. 
Therefore, a good way to calculate the intermodulation noise in practical 
transmission systems and to analyze effects of the transmission characteris- 
tics in trunk on the FM multiplex telephony signal of radio relay systems are 
provided. (Paper received July 1984, finalized February 1985.) 
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INPUT METHODS OF CHINESE INFORMATION PROCESSING SYSTEM 


Beijing DIANZI XUEBAO [ACTA ELECTRONICA SINICA] in Chinese Vol 14 No l, 
Jan 86 pp 99-107 


[English abstract of article by Zhao Pozhang [6392 3789 3864] of the Beijing 
Institute of Informationand Control] 


[Text] The five-layer structure model for Chinese information processing 
system is proposed, andthe features of its input method are discussed. Also 
presented are the idea of general purpose Chinese character input coding 
module and the technical feature of natural input method--character recogni- 
tion and speech recognition briefly. (Paper received May 1985, finalized 
September 1985.) 
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OBSERVATION ON RELATIONSHIP BETWEEN MEASUREMENT OF BLOOD FIBRONECTIN AND 
CLINICAL COURSE OF BURNS 


Beijing ZHONGHUA WAIKE ZAZHI [CHINESE JOURNAL OF SURGERY] in Chinese Vol 24 
No 4, 22 Apr 86 pp 222-224, 253 


[English abstract of article by Huang Wenhua [7806 2429 5478], Huang Xingyu 
[7806 5281 3768], and Li Ao [7821 7663] of the First Affiliated Hospital of 
the Third Military Medical College] 


[Text] Blood fibronectin (FN) is a nonspecific opsonic glycoprotein. In the 
phagocytic and clearing function of reticuloendothelial system, it acts as an 
important component. It is one of the useful indexes for monitoring critical 
patients suffering from severe trauma, shock, or sepsis. This paper presents 
the kinetic changes of blood FN, measured by agglutination assay, in 32 cases 
of burned patients. The preliminary impression was that the changes of FN 
titre closely correlated with the conditions of patients. In early postburn 
period (1 st wk), blood FN reduced markedly, its level paralleled with the 
severity of injury. In late postburn period, drop of blood FN titre to low 
level indicated the development of severe systemic and/or local sepsis as 
well es other complications. 
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